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I.  Numerical-analytical  model  of  penetration  of  long  elastically  deformable  projectiles  into 
semi-infinite  targets 


3 


Introduction 

Penetration  of  non-deformable  projectiles  in  continuum  with  various  rheological  properties  has  been  of 
interest  to  researchers  for  a  long  time.  The  first  modeling  representations  of  penetration  were  formulated  in 
XVIII-XIX  centuries  in  Euler's,  Poncelef  s,  Wuich’s  works,  etc.  Analysis  of  these  results  can  be  found  in  A.J. 
Sagomonjan's  monographs  [1,2]. 

In  conjunction  with  development  of  more  exact  and  effective  technical  means  in  last  two  decades  the 
interest  to  this  problem  has  considerably  increased,  and  this  is  proved  by  works  of  Voejkova  and  Sagomonjan 
(1985)  [3],  Alojan  (1985)  [4],  Liapykhin  et  al.  (1993)  [5],  Bahrah  et  al.  (1992)  [6],  Forrestal  et  al.  (1988)  [7], 
(1992)  [8],  (2000)  [9],  Dikshit  and  Sundararajan  (1992)  [10],  Piekutowski  et  al.  (1999)  [11],  Warren  and 
Forrestal  (1998)  [12],  Warren  (2000)  [13],  Yossifon  et  al.  (2001)  [14],  Chen  and  Fi  (2002)  [15]. 

The  analysis  of  modern  state  of  the  problem  of  analytical  modeling  of  high-velocity  penetration  of  non- 
deformable  projectiles  in  target  can  be  found  in  works  by  Forrestal  et  al.  [7-9],  Warren  and  Forrestal  [12],  Yarin 
et  al.  [16,17],  Yossifon  [14],  Chen  and  Fi  [15].  From  this  analysis  follows  that  at  present  there  is  a  deficiency  of 
relatively  simple  analytical  models  using  natural  physical  and  geometrical  parameters  of  projectile  and  target 
and  with  a  small  number  of  fitting  parameters. 

In  this  study  we  have  built  and  investigated  a  new  model  of  penetration  of  non-deformable  projectiles  of 
various  shapes  in  elastic-plastic  and  elastic-brittle  materials. 

The  closest  to  the  proposed  model  are  models  by  Yarin  et  al.  [16],  Roisman  et  al.  [17],  Yossifon  et  al. 
[14]  which  use  the  method  of  peculiarities  where  equations  of  motion  in  elastic-plastic  medium  have  a  form  of 
the  Cauchy  problem: 

(w -r  -r  C(x,i),  x(0)  =  0,  i(0)  =  -f/o.  (1) 

Here  Uo>0  -  initial  projectile  velocity,  m  -  its  mass,  x(t)  -  location  of  the  tip  of  the  projectile,  and 
functions  A,  B,  C  have  natural  mechanical  meaning  {A  -  target  added  mass,  B  -  drag  coefficient,  C  -  drag 
associated  with  plasticity).  These  functions  are  found  analytically  on  the  base  of  approximate  solution  of  motion 
equations  of  elastic-plastic  medium  with  appropriate  boundary  conditions.  Equations  (1)  use  exact  velocity 
potential  of  ideal  fluid  in  infinite  domain  for  a  projectile  which  shape  is  an  ovoid  of  Rankine. 

The  comparison  of  full  penetration  depth  of  long  non-deformable  projectile,  determined  by  our  model 
for  different  velocities  of  impact  with  known  experimental  data  shows  that  the  experimental  values  of 
penetration  depth  appear  less  than  theoretical  ones.  Similar  results  are  demonstrated  by  the  model  suggested  in 
[7,8]  and  based  on  a  use  of  velocity  field  in  material  at  expansion  of  spherical  cavity.  To  adjust  the  theoretical 
results  to  experimental  ones  the  authors  in  these  works  use  a  friction  coefficient.  Eater,  on  the  basis  of  the 
analysis  of  FEM-based  calculations  and  experimental  data  [12],  they  neglect  the  influence  of  friction  on 
penetration  depth  and  develop  the  model  that  accounts  for  deformation  rate  in  rheological  equations  of  target 
material.  Analyzing  the  influence  of  hardness  of  projectile  material  on  experimental  values  of  penetration  depth 
in  work  [9]  it  is  noted  that  harder  projectiles  have  a  slightly  higher  penetration  depth  than  less  hard  ones. 
However,  all  experimental  values  of  penetration  depths  are  lower  than  theoretical  depths  for  hard  projectiles. 

These  facts  led  us  to  the  necessity  to  take  into  account  other  possible  effects,  in  particular,  elastic 
deformation  of  projectile.  In  present  study  we  investigate  in  uniaxial  approximation  the  influence  of  forced 
elastic  oscillations  of  projectile  on  key  parameters  of  penetration  process  (penetration  depth  and  velocity).  It  is 
shown,  that  elastic  oscillations  of  projectile  (usually  neglected)  have  noticeable  influence  on  penetration. 

We  also  evaluate  a  critical  velocity  of  impact  at  which  the  projectile  undergoes  plastic  deformation 
causing  irreversible  changes  of  its  shape  and  instability  of  its  movement.  The  investigation  of  rigid-body  motion 
and  further  analysis  of  forced  elastic  oscillations  deals  with  impact  velocities  less  then  this  critical  velocity. 

The  main  purpose  of  the  present  modeling  investigation  is: 

-  to  build  motion  equations  on  the  basis  of  velocity  field  in  the  target  determined  by  a  real  projectile 

shape; 

-  to  obtain  deeper  insight  into  the  nature  of  target  resistance  to  penetration  of  non-deformable  and 
elastically  deformable  projectiles; 

-  to  predict  the  dependence  of  projectile  velocity  vs.  given  geometrical  and  material  parameters  of  projectile 
and  target; 

-  to  investigate  the  influence  of  projectile  shape  and  its  elastic  oscillations  on  penetration; 

-  to  compare  results  of  model  predictions  with  known  experimental  and  calculated  results. 
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1.  Basic  model  equations  of  penetration  of  non-deformable  projectiles  into  semi-infinite  target 


1.1.  Geometrical  scheme  of  penetration  model  of  non-deformable  projectile 

Let’s  consider  a  normal  penetration  of  axisymmetric  projectile  in  semi-infinite  target  (see  Fig.l). 

In  Fig.  1  we  accept  the  following  notations  and  assumptions:  Oi  —  lateral  surface  of  projectile  (i.  e. 
projectile  surface  without  its  rear  surface),  02  =  {x  =  0,  r  >  dl2)  —  target  free  surface  (assumed  flat),  03  =  {0  <  x 
<  (P  -  L),  r  =  d/2}  —  free  surface  of  entry  hole  (assumed  cylindrical).  The  direction  of  projectile  motion 
coincides  with  axis  x  which  is  its  axis  of  symmetry.  AtP  <L  surface  03  is  absent. 

t 

Main  model  parameters  are:  u{t)  —  projectile  velocity;  L  —  projectile  length;  P{t)-^u  {x)dx  — 

0 


penetration  depth;  d  —  projectile  diameter;  Us  —  striking  (impact)  velocity. 


Fig.  1  Scheme  of  normal  penetration  of  axisymmetric  elongated  non-deformable  projectile  in  cylindrical  system 

of  coordinates  OrOx 


1.2.  Velocity  field  in  target 

Velocity  field  in  the  vicinity  of  projectile  is  approximated  by  the  field  of  velocities  v  for  irrotational 
motion  of  ideal  incompressible  fluid.  Defining  v  we  assume,  that  a  free  surface  of  target  02  =  {x  =  0,  r  >  dl2} 
remains  fiat  during  penetration,  and  a  free  surface  of  entry  hole  03  formed  by  a  tail  part  of  projectile  is 
cylindrical:  03  =  {r  =  dl2,  0<x<{P-  L)},  (see  Fig.  1). 

As  it  is  known  [18,19]  the  equations  of  non-stationary  irrotational  movement  of  ideal  incompressible 
fluid  with  density  p  in  case  of  potential  flow  have  the  following  form: 

V^cJ)  =  0,  Laplace  equation  (2) 

p{d<^/dt)  +  p/2|grad(I)|^  -I- p  Lagrange-Cauchy  integral  (3) 

where  <P(x,r,t,P)  -  velocity  potential  (i.e.  v  =  -gradO),  p  -  static  pressure  component, y(0  -  arbitrary  function  of 
time,  independent  from  coordinates  and  identical  over  the  whole  volume  of  fluid.  Function  y(0  is  found  from 
boundary  conditions. 

First  equation  (2)  is  a  condition  of  incompressibility  of  fluid  (or  the  continuity  equation),  second 
equation  (3)  -  the  first  integral  of  Euler’s  equation  of  motion  which  is  known  as  Lagrange-Cauchy  integral 
(sometimes  it  is  called  Bemoulli-Cauchy  integral).  Equation  (3)  is  conservation  of  momentum  and 
grad[p(5(I)/5f)  +  p/2|grad(I)p  +  p}  =  diva'  =  grad(/(f))  where  a'  is  a  deviatoric  part  of  the  Cauchy  stress  as 
described  by  Yarin  et  al  [16].  If  Eagrange-Cauchy  integral  (3)  is  considered  on  projectile  surface  left  part  of 
equation  (3)  defines  a  normal  pressure  pc  of  target  material  on  projectile  surface,  i.Q.pc  =  p{d(^ldf)  +  p/2|grad(I)|^ 

+p. 

Due  to  the  assumptions  made  and  because  equation  (2)  does  not  contain  partial  derivative  with  respect 
to  time,  time  t  enters  into  (2)  as  a  parameter  which  influence  is  manifested  through  the  boundary  conditions 
containing  geometrical  characteristics  of  projectile  and  its  velocity  u(t).  Therefore,  both  time  and  depth  of 
penetration  P(t)  are  listed  in  the  arguments  of  potential  (J). 

Thus,  for  potential  (J)  we  have  the  following  Neumann  boundary  conditions: 

d<^/dn  =  -  u{t)nx  at  M(r,x)  eai,  d<^/dn  =  0  at  M(r,x)  e(a2  u  a3), 

(I>=  0(l/i?)— >0  when  R  =  +x^  — >•  00  , 


(4) 
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where  ai,  a2,  as  —  surfaces  as  indicated  earlier,  —  cosine  of  angle  a  between  axis  x  and  normal  to  projectile 
surface  ai,  (see  Fig.  1).  Condition  (4)  is  the  condition  of  impenetrability  of  target  material  across  projectile 
surface  ai,  cylindrical  surface  as,  and  target  surface  as. 

From  the  linearity  of  boundary  problem  of  Neumann  (2),  (4)  follows,  that  velocity  potential 
may  be  presented  as 

<^{x,r,t,P)  =  u{t)(^{x,r,P).  (5) 

Here  function  (p(x,r,F’)  is  a  solution  of  Laplace  equation 

V^cp  =  0  (6) 

with  boundary  conditions 

5(p/5n  =  -«;£  when M(x,r)e ai,  0(p/5n=O  atM(x,r)e(a2  u  as),  (p  =  0(l/f?)— >0  whenf?  — >oo.  (7) 

Let's  notice,  that  (p(x,  r,  P)  depends  on  P  as  on  parameter,  which,  in  its  turn,  depends  on  time  t. 

Thus,  the  field  of  velocities  v  in  the  proposed  model  is  defined  by  an  expression 

V  =  -u{t)  grad  (p(x,  r,  P).  (8) 

Obviously,  velocity  field  (8)  is  approximate,  since  during  penetration  a  free  surface  of  target  as  does  not  remain 
plane  and  a  free  surface  of  input  hole  as  is  not  a  circular  cylinder  (see  Fig.  1).  In  a  complete  formulation  of 
penetration  problem  of  rigid  projectile  the  shape  of  these  surfaces  has  yet  to  be  determined. 

However,  according  to  the  experiments  (for  example,  Piekutowski  et  al,  1999  [11])  the  deviations  from 
surfaces  as  and  as  accepted  in  the  proposed  model  are  insignificant  (in  comparison  with  a  diameter  of 
projectile),  especially  it  concerns  plastic  materials  (see,  for  example.  X-ray  photo  from  the  mentioned  work.  Fig. 
1).  Thus,  the  accepted  assumptions  as  to  surfaces  as  and  as,  in  our  point  of  view,  insignificantly  distort  actual 
field  of  velocities  v ,  but  considerably  simplify  the  investigation  of  this  field  and  emphasize  its  main  features. 

Let's  note,  that  in  the  proposed  model  the  function  (I)  =  wcp  is  a  velocity  potential  of  ideal  incompressible 
fluid  for  which  there  are  no  friction  forces  (tangential  stresses)  between  target  material  and  projectile  surface. 
This  approximation,  which  neglects  tangential  stresses  on  projectile  surface,  is  a  common  one  in  the  majority  of 
computing  schemes.  As  was  mentioned  earlier,  in  [12]  on  the  base  of  computational  modeling  and  experimental 
evidence  authors  rejected  the  use  of  sliding  friction. 

Because  of  this,  we  believe  that  a  prediction  of  integral  values  of  penetration  process  (such  as  a  force 
that  decelerates  projectile,  penetration  depth  and  velocity)  which  use  velocity  field  (8)  can  be  quite  satisfactory. 
Moreover,  as  the  main  influence  on  projectile  motion  is  rendered  by  potential  (I)  =  wcp  and  velocity  v  =  -wgradcp 
on  the  surface  of  a  head  part  of  projectile,  then  the  influence  of  boundary  conditions  (7)  on  as  and  as  would 
decrease  as  penetration  depth  P  increases  and  become  insignificant  when  P  would  exceed  several  projectile 
head  lengths. 

1.3.  Ballistic  integro  differential  equations  of  projectile  motion 

As  in  [20],  for  building  of  projectile  motion  equation  we  use  Lagrange-Cauchy  (3).  If  we  substitute 

velocity  potential  (5)  (O  =  M(t)(p(M,T’),  P-\u{x)dx,  =  cp— -t ,  gradO  =  wgradcp)  in  Lagrange- 

“  dt  dt 

Cauchy  integral  (3)  and  consider  this  integral  on  projectile  surface,  then  a  left-hand  part  of  (3)  would  define  a 
contact  pressure  pc(M,P)  of  target  material  on  projectile.  Here  coefficient  cp),  takes  into  account  change  of  the 
target  boundary  with  penetration,  i.  e.  domain  of  definition  of  velocity  field  (8).  Therein  lies  its  mechanical 
significance.  Further  we  assume  that  in  (3)  static  components  of  pressure  on  projectile  surface  p  =  Rt  where  Rt 
—  pressure  in  a  spherical  cavity  during  its  quasi-static  expansion,  i.e.  here  we  use  the  analytical  solution  of  a 
problem  of  quasi-static  expansion  of  spherical  cavity  in  target  material  which  takes  into  account  elastic  and 
plastic  (or  brittle)  response  of  material.  Rt  is  constant  and  depends  on  mechanical  properties  of  target  material 
only  (see  (16-19)).  The  use  of  such  solution  is  caused  by  a  fact  that  analytical  solution  of  quasi-static 
indentation  of  rigid  projectile  of  real  shape  in  target  material  generally  is  not  possible.  Besides,  since  target 
material  at  this  indentation  flows  or  fractures  and  relaxation  of  stresses  takes  place,  value  Rt  (as  it  will  be  shown 
further)  enters  into  the  equation  of  projectile  motion  as  integrand.  This,  in  our  opinion,  makes  the  adopted 
approximations  justified  and  sufficient. 

Thus,  on  the  basis  of  Lagrange-Cauchy  integral  (3)  for  normal  dynamic  pressure  of  target  material 

/  \  50  1  I  i2 

p^iM,P)  =  p - -t  —  p  gradO  +  p  on  projectile  in  point  Mon  its  surface  we  obtain 

dt  2 

2 

p,  {M,P)  =  p(p(M,P)^  +  ||grad(p(M,P)|"  +  2(p>  (M,P)  +  R^ , 


(9) 
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where  P  -^u{x^dx  — projectile  penetration  depth. 


A  full  decelerating  force  N  acting  on  projectile  is  obtained  by  integration  on  a  contact  surface  of 
projectile  and  target  a{P)  of  a  projection  (-pc(M,P)-nJ^M,P))  of  contact  pressure  (9)  onto  axis  x.  Here  nJ^M,P)  — 
a  cosine  of  an  angle  which  is  formed  by  unit  normal  to  projectile  surface  in  point  M  with  axis  x  (see  Fig.  1). 

dti 

Equating  a  complete  force  of  resistance  to  projectile  A  to  its  inertia  force  m —  (m  —  projectile  mass) 

dt 

the  ballistic  equation  of  motion  of  non-deformable  projectile  is  obtained: 

nil 

m  —  =N,  M(0)  =  f/o,  (10) 

dt 

where  Uo=  —  impact  velocity. 

Accounting  for  (9),  equation  ( 1 0)  can  be  written  as 


m^  =  -pA{P)^-B{P)^-C{P)R,,u{0)  =  Uo.  (11) 

ai’  HL  ^ 

In  (1 1)  the  coefficients  A,  B,  C  depend  on  projectile  shape,  penetration  depth  P  and  are  defined  by  the  formulas 
A(E)=  j  i^{M,P)n,{M,P)daM  ,  C{P)^  j  (12) 

a(P)  a(P) 


B{P)^  j  |grad(p(M,E)|  +  2(p), (M,E)  . 


In  these  formulas  a  contact  surface  a(E’)  of  projectile  and  target  in  a  moment  of  time  t  is  approximately  defined 
by  the  inequality 


a(E)  =  ^M:p(p(M,E)— +  ^|grad(p(M,E)|  +2(pp(M,E)j^  +  i?,  >0  . 

If  projectile  surface  is  formed  by  a  rotation  of  line  r  =  r(x)  around  axis  x  (see  Fig.  2),  then  in  (12) 

=  ^  ,  dotf  -2TirJl  +  (r')^dx, 


and  these  formulas  can  be  written  as: 


A(P)  -  -2tc|  (p(x,r(x),E’)r(x)r'(x)G?x,  ~  ri^x)r' i^x)dx. 


(E)  = -27r|Mgrad(p(x,r(x),E’)  +  2(pp (x,r(x),E’)  r(x)r'(x)G&c, 


where  a  stroke  denotes  derivative  with  respect  to  x,  and  Ip  —  intervals  on  axis  x  which  define  a  region  of  contact 
of  projectile  and  target  material: 

/p  =  |x:p(p(x,r(x),E’)^  +  |^|grad(p(x,r(x),E’)|^  +2(pp(x,r(x),E’)j^^  +  f?,  >o| .  (15) 

Fig.  2  illustrates  one  contact  interval  lp  =  [Xc,P]. 


Fig.  2  To  the  definition  of  contact  surface  between  projectile  and  target  material;  C  —  point  of  separation  of 

target  material  from  projectile  surface 
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Inequalities  (13)  and  (15),  approximately  determining  contact  surface  a(/’)  between  projectile  and  target 
material,  are  more  detailed  form  of  the  condition  pc{M)  >  0.  Regions  of  separation  of  target  material  from 
projectile  surface  a(R)  are  approximately  defined  by  a  condition:  a(R)  =  {M:  pc<0}. 

As  it  was  noted  earlier,  the  value  of  static  penetration  resistance  of  target  material  Rt  is  defined  by  the 
analytical  solution  of  the  problem  of  expansion  of  a  spherical  cavity.  At  present  there  exist  a  number  of 
solutions  for  both  brittle  and  plastic  materials.  Several  of  them  that  are  used  further  we  list  below. 

For  elastic-ideally  plastic  and  elastic -brittle  materials  the  following  formulas  can  be  used  respectively 

[21-23]: 


R,-h 

t  3 


1-tln 


V 


37(l-v) 


Rt-Y 


-i2a/3 


37(l-(l-v)^a//(27)) 


(16) 


(17) 


Here  E,  v  are  Young’s  modulus  and  Poisson’s  ratio,  Y —  compressive  yield  strength  for  (16)  and  strength  at 
uniaxial  compression  for  (17),  a/ —  tensile  strength,  a  —  characteristic  of  fractured  media  [23].  To  account  for 
a  post-yield  strain-hardening  of  material  the  following  formulas  can  be  used  [8,12] 


t  3 


r 


1-t 


V 


2E 

37 


/ 


l-'iYl(2E) 


(-Inx)^ 
1  -X 


dx , 


(18) 


where  n  —  strain  hardening  exponent. 

For  penetration  of  in  limestone  Frew  [24]  established  that  penetration  resistance  depends  on  penetrator 
dimensions: 

= 'F -t  (|)  (flo/a),  (19) 

where  'F  =  607  MPa,  (|)  =  86  MPa,  2  ao  =  25.4  mm,  a  —  projectile  radius  (mm). 

1.4.  Determination  of  velocity  potential  cp  by  method  of  boundary  integral  equations  [25] 

The  solution  of  the  Neumann  problem  (6),  (7)  is  sought  as  a  potential  of  a  simple  layer 

^dS{M,),M^[DKjS),  M,gS,  (20) 

in  which  S  =  (5'u5")  —  a  surface  symmetric  to  plane x  =  0,  and  D  —  infinite  area  with  a  boundary  S  ; 
R{M,M\)  —  distance  between  a  fixed  point  of  observation  M(x,r,0)  and  a  source  point  of  integration  Mi(xi,ri,0i) 
(see  Fig.  3);  p(Mi)  —  the  unknown  density  of  simple  sources  distributed  on  surface  S  . 

The  density  p(Mi)  is  a  solution  of  the  following  boundary  integral  equation 
1  I-  cos( 


i{R,n\  1 

R  (M,Mi)  2ti 


(21) 


g(M)  = 


MeajUaj, 

.X 


(22  a-b) 


[O,  MGa3Ua3, 

where  vector  R  is  directed  from  point  M  to  Mi,  — unit  normal  to  S  in  point  M  (internal  with  respect  to 

D ),  nx{M)  —  a  projection  of  a  normal  n  onto  axis  x. 


D’ 


S'^  (a'lUo'j) 


L 


Plane  of  symmetry 

(x  =  0) 


Fig.  3  Geometrical  scheme  of  the  domain  of  definition  D-DkjD'  for  potential  cp.  Region  D'  and  surface 
S'  =  (a'iUa'3)  are  symmetric  to  region  D  and  surface  S  =  (01^03)  with  respect  to  plane  x  =  0. 


If  surface  S  is  formed  by  a  rotation  of  line  r  =  r(x),  x&[-P,P]  around  axis  x  and  symmetry  to  a  plane  x 
=  0  and  axis  x  is  taken  into  account  as  well,  then  boundary  equation  (21)  can  be  written  as 

p(x)  +  —  [\  Ki^X,Xy)  +  ki^X,-Xij\\li^Xx)dXi  =  \  X,  Xi  G  [0,R], 

TT  j  1-  -I  2tc 


g(x)- 


0,  xg[0,R-L[, 

xg[0,R]\[0,R-L[. 

[  Vl  +  r'2 


Here  we  have  the  following  notations: 

^£'(m)[(x-xi)n^  +(r-ri)n^]  E{m)-  K{m)i 


^(x,xi)  = 


(a  -b)\la  +  b 


2r\ja  +  b 


-  +(^x  -  XiY ,  b  =  2rr\,  r  =  r(x),  rj  =  r(xi),  m  - 


2b 
a  +  b 


,  K  —  full  elliptic  integral  of  the  first  order,  E  — 


full  elliptic  integral  of  the  second  order  [26],  nr{M)  —  a  projection  of  normal  n  to  axis  r  (see.  Fig.  3). 


For  components  of  velocity  =  -uit) — ,  =  -u(t) —  in  directions  x,  r  (see  Fig.  1),  respectively,  we 

dx  dr 


obtain 


,{x,r)  =  4u{t)  j  f  fi 

2p[a-  b)\ia  +  b 


fXi 


v^{x,r)  =  4u{t)  j 


-p 


2r(a  -b)\ja  +  b  2r^a  +  b 


rj-Jl  +  r/^p(xi)(ixi  .  (23) 


2.  An  approximate  account  of  elastic  oscillations  of  long  projectiles  and  estimation  of  value  of 
impact  velocity  at  which  a  projectile  undergoes  a  plastic  flow 

Further,  while  evaluating  the  influence  of  elastic  oscillations  of  long  projectile  on  its  penetration,  we 
neglect  local  elastic  deformations  of  projectile  nose  and  assume  that  the  projectile  has  the  shape  of  circular 
cylinder  with  a  diameter  d  and  length  L.  These  oscillations  are  considered  separately  from  projectile  motion  as  a 
rigid  body  and  the  estimate  of  their  influence  on  the  transportation  velocity  of  penetrator  is  made  via  energy 
loss.  To  evaluate  the  influence  of  forced  elastic  oscillations  of  projectile  we  use  a  solution  of  a  problem  of 
longitudinal  oscillations  of  a  rod  [27,28]  with  one  free  end  (x  =  0),  and  compressive  stresses  (-ao)  instantly 
applied  to  another  (x  =  L)  (at  the  moment  of  time  f  =  0): 
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I  2ao-^ (-1)^  r,  incnt  iTix\  . 

W  =  -  +  1-cos^  cos—  ,  i=l,2,3,...,  (24) 

[2ppL  I  L  }  Lj 

where  w{x,t)  —  relative  displacement  of  cross-section  of  a  rod  x  at  the  moment  of  time  t,  —  density  of  rod 
material,  Ep  —  Young’s  modulus  of  rod  material,  Cq  =  ^E^jpp  —  sound  speed  in  projectile.  The  influence  of 

cross-section  oscillations  is  neglected.  We  believe  the  use  of  this  solution  is  justified  since  (as  the  calculations 
show,  see  Fig.8)  a  penetration  of  long  rigid  projectile  occurs  at  almost  constant  acceleration  [dujdt)  and,  thus, 

d 

at  almost  constant  braking  force  N  --m —  =  — - >0,  acting  on  projectile  (here  u  —  solution  of  initial 

dt  4 

boundary  problem  (11)).  Besides,  the  value  of  force  N  weakly  changes  at  small  variations  of  initial  condition 
m(0)  =  Uo  in  (11). 

Let's  note,  that  in  this  paragraph  for  convenience  it  is  supposed,  that  the  force  compressing  projectile  N 
is  positive,  unlike  (10). 

From  (24)  for  displacement  w  and  velocity  w  atx  =  L  (i.e.  for  w  and  w  at  projectile  tip  caused  by  its 
forced  oscillations  with  a  period  x  =  2L/co  )  we  obtain 

—  (25) 

Ep  dt  Ep 

4N 

Here  Gq  =  — y  and  braking  force  >  0  is  a  little  changing  value  {N  »  const  through  a  whole  depth  of 
nd^ 

penetration);  vectors  w  and  w  are  directed  opposite  to  the  direction  of  projectile  motion.  Therefore  the  initial 
velocity  of  projectile  penetration  Uo  in  the  target  is 

(26) 

Ep  nd  Ep 

where  Us  —  projectile  velocity  before  its  meeting  with  target,  i.e.  impact  velocity.  For  a  non-deformable 
projectile  Uq  =  Us. 

On  the  basis  of  (24)  a  longitudinal  deformation  of  projectile  dwjdx  (caused  by  oscillations)  at  the 
moment  of  time  t  is  defined  by 


dw  _  2gq  (  l)  j 
dx  nEp  ^  i 


mcU  1  .  inx 


sin - ,  /=  1,2,3, 

L 


For  the  moments  of  time  t-—{2i-\'\  this  deformation  —  = --^=  const  <  0  at  all  values  xe[0,L].  Therefore 

Co  dx  Ep 

during  these  moments  of  time  projectile  is  entirely  compressed  by  stresses  (-ao),  i.e.  in  each  cross-section  of 
projectile  x  there  exists  compressing  stress  (-ao). 

1  2 

During  an  impact  the  energy  of  projectile  consists  from:  a  kinetic  energy  T  -—mu  (u  —  projectile 

transportation  velocity);  elastic  energy  Wi  caused  by  lateral  compression  of  projectile;  energy  of  elastic 
oscillations  of  projectile  W2,  caused  by  sudden  application  of  braking  force  N. 

The  elastic  energy  Wi  associated  with  lateral  transversal  compression  of  projectile  by  pressure 


Pi  -  p^iyM  ,P^sma-  p(p(M,T’)— +  |^|grad(p(M,T’)|  +2(pp(M,T’)j^— +  sina, 

where  point  M  is  on  a  lateral  surface  of  projectile  and  a  —  an  angle  which  forms  a  normal  to  projectile  surface 
in  point  M  with  its  axis,  pc  —  normal  pressure  of  target  material  on  projectile  surface,  see  (9),  is  approximately 
estimated  by  the  formula 

[  e{pi)dx,  P<L, 

(28) 

j^e(pi)dx,  P>L. 
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nd^\\-v  p\ 

Here  e{^pi)  = - pi  [M,P^  —  energy  of  deformation  of  projectile  cross-section  in  point  M  at 

4Ep 

penetration  depth  P,  determined  in  the  assumption,  that  normal  stress  in  projectile  cross-sections  is  absent  (v^  — 
Poisson  ratio  for  projectile  material).  Estimating pi  we  assume  sina  =  1. 

Energy  of  forced  elastic  longitudinal  oscillations  of  projectile  W2  accumulated  to  the  moment  of  time  t 
(work  of  drag  force  spent  on  elastic  oscillations)  on  the  basis  of  (25)  is  estimated  by  value 

=  =  (29) 

2  nd  Ep  ^Ep 

which  coincide  with  value  [27]  where  it  is  represented  by  a  convergent  series.  From  (29)  can  be  seen  that  the 
energy  of  forced  oscillations  increases  proportionally  to  time  while  braking  force  remains  constant. 

AN 

Since  p;  <  Gq  =  — y ,  then  for  energy  W\  we  obtain  an  estimation 


nd 


< 


4(\-Vp)N^L 


nd^E, 


From  this  and  from  (29)  at  cot  »  2L  we  obtain  W\  «  W2  and  we  can  neglect  the  energy  of  transverse 
compression.  Further  calculations  confirm  that  we  can  neglect  this  energy. 

Solution  of  equation  (1 1)  at  initial  condition  (26)  we  denote  as  u  and  further  refer  to  it  as  “a  velocity  of 
penetration  of  quasi-rigid  projectile”.  This  can  be  justified  by  the  fact  that  initial  condition  (26)  for  the  equation 
of  motion  of  non-deformable  projectile  (1 1)  takes  into  account  projectile  elastic  deformations. 

1  _ 2 

Thus,  between  kinetic  energy  of  quasi-rigid  projectile  T  -—mu  and  complete  energy  of  deformable 
1  2 

projectile  —mu  +W1  +  W2  there  is  a  relationship 

1  _ 2  1  2 

—mu  -—mu  +IT1+IT7, 

2  2  1 

which  defines  penetration  velocity  u  of  deformable  projectile  through  quasi-rigid  projectile  velocity  u 


u  =  ujl- 


2{W,  +  W2) 


-2 

mu 


(30) 


du 


where  values  Wi,  i  =  1,2,  are  calculated  by  formulas  (28),  (29)  and  depend  on  braking  force  N  -  -m — ,  which 

dt 

is  approximately  defined  by  Cauchy  problem  (11)  and  (as  it  was  discussed  earlier)  weakly  changes  during  an 

du 

impact.  Thus,  with  a  sufficient  accuracy  one  may  also  assume  N  »  -m  —  .  If  we  neglect  energy  ,  then 

dt 

formula  (30)  takes  simple  form 


u-ujl- 


2W2 

-2 

mu 


(30') 


From  (30)  follows,  that  a  motion  of  projectile  (as  a  solid  body)  with  transportation  velocity  u  takes  place 
until  the  moment  of  time  t  when  the  following  equality  comes  true 

2{Wi+W2)  =  mu^ . 

1  2 

It  means  that  a  penetration  of  projectile  stops  when  its  kinetic  energy  —mu  becomes  zero  and  it  has  only  the 

energy  of  its  elastic  deformation.  After  this  moment  a  transportation  motion  of  projectile  stops  and  it  performs 
damped  oscillations. 

Having  defined  a  transportation  motion  of  projectile  u(t)  by  (30)  for  a  complete  actual  penetration  depth 
P  we  obtain 


P  =  ^u{xyix , 


and  for  a  correction  AP  to  a  final  penetration  depth  of  quasi-rigid  projectile  P  we  have 
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AP^P-P^ 


I  u  [x)dx  -  I  M  {x)dx  <  0  , 


where  t  —  time  of  impact  of  quasi-rigid  projectile  (t  <t  ). 

Let  evaluate  a  threshold  value  of  impact  velocity  Us  =  Usc  at  which  a  projectile  undergoes  plastic 
deformations  causing  irreversible  changes  of  its  shape  and  instability  of  motion. 

Since  ao  —  peak  value  of  an  axial  stress  in  projectile  (if  a  projectile  does  not  lose  stability  at 
compression)  the  plastic  deformations  in  it  start  to  develop  under  a  condition 

4N 

(31) 

na 

where  Yp  —  dynamic  yield  stress  of  projectile  material. 

,  f  2  '\ 

Cltl  fft  Otl 

If  we  take  into  account  N  -  -m  —  = - ; — -  B(P) - +  C(P)Rf  (see  equation  (11)),  then  from 

dt  m  +  pA{P\  ^  ’  2  V  I  V 

(3 1)  we  obtain  the  following  equation  for  projectile  velocity  u  at  which  the  plastic  deformations  appear  in  it 


ffl  +  -y  ,  (32) 

where  A{P),  B{P),  C{P)  —  coefficients  of  the  equation  (11)  and  m  =  pA(Pj  —  an  added  mass.  Let’s  note  here, 

that  the  coefficients  A,  B,  C  depend  on  penetration  depth  P  (as  parameter)  and  shape  of  the  projectile  nose.  At 
sufficiently  large  values  of  P  (for  projectiles  discussed  further  P  >  L^,  Lf,  —  length  of  a  head  part  of  projectile) 
they  are  constant  (do  not  depend  on  P  and  depend  only  on  projectile  shape).  Therefore  for  these  values  P  and 
the  right-hand  part  of  equality  (32)  is  a  constant  value.  Besides,  as  further  calculations  show,  this  constancy  is 

2L 

almost  established  by  the  moment  of  returning  of  a  front  of  elastic  longitudinal  wave  (t  -  — )  in  projectile  to  a 

Co 

contact  surface.  Therefore  further  factors  A,  B,  C  are  considered  as  constants,  and  argument  P  is  not  referred  to 
in  formulas. 

As  the  largest  value  of  velocity  u  is  observed  at  the  beginning  of  impact,  when 


(see  (26)  and  (31)),  for  critical  velocity  of  projectile  Usc  till  its  meeting  with  a  target 


(i.e.  for  critical  velocity  of  impact)  from  (32)  we  obtain  the  estimate 

Yp<^o  .  \nd^  f r,  .  T 


2p5^^  mj  ^  nd^ 


and  projectile  flows  under  a  condition  Us  >  Usc-  Thus,  plastic  deformations  in  projectile  develop  from  the 
beginning  of  impact.  In  formula  (33)  radicand  has  to  be  positive.  If  this  expression  is  negative,  i.e.  the  opposite 
inequality  is  true 


i.%-TL«,so 

mJ  ^  Tzd 


which  means  that  the  target  is  of  the  same  strength  as  the  projectile  or  even  stronger,  then  critical  velocity  of 
impact  is 

u 


and  projectile  flows  under  a  condition  of  Us  >  Usc  which  coincides  with  known  condition  [29].  Thus,  the 
condition  Us  >  Usc  where  Usc  is  defined  by  (33),  is  a  generalization  of  this  known  condition  of  transition  of  the 
projectile  material  into  a  plastic  state. 

0^4  _  4  _  4  C 

For  long  projectiles  when  —  = - « 1 ,  A-  — -  and  C  =  — -  » 1 ,  inequality  (34)  simplifies  and  has 

m  PpL  Tid^  nd^ 

a  simple  mechanical  meaning:  R,  >  Yp. 
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If  we  introduce  “reduced  strengths”  of  projectile  and  target,  respectively, 

pA 


Rf  —  CRf , 


1  +  - 


9pL  j 


(35) 


-  AB 

and  denote  B  - 


nd^ 


,  then  formula  (33)  for  critical  velocity  can  be  rewritten  as 


YpCo 


YpCo 


■+  -^{y-r,] 
]Ipb''  ^  ' 


Yp<Rt, 


(36) 


and  is  analyzable  by  elementary  means.  Indeed:  if  value  Yp  increases,  then  critical  velocity  Usc  increases  also;  if 
strength  of  target  Ri  falls,  this  speed  increases;  if  sharpness  of  projectile  nose  increases  the  coefficient  B 
decreases  (for  spherical  (vj/  =  0.5)  nose  B  =  0.58,  for  ogival  (\|/  =  3)  nose  5  =  0.1 1)  and,  hence,  Usc  increases;  if 
bluntness  of  the  nose  increases  the  coefficient  A  increases  (for  ogival  nose  A  =  0.0010,  for  spherical  A  = 
0.0015)  and,  therefore,  the  reduced  strength  of  projectile  Yp  (associated  with  the  added  mass  m-pA,  see  (35)) 

increases.  Formula  (36)  can  be  used  to  define  a  dynamic  yield  stress  of  projectile  material  Yp  if  other  values 
entering  into  it  are  known.  See  example  in  appendix. 

Note,  that  the  case  of  stability  loss  of  long  projectile  can  be  also  considered  under  the  proposed  scheme 
of  investigation  of  plastic  deformations  of  projectile.  Indeed,  in  this  case  in  the  main  formulas  (31)  —  (36)  it  is 
necessary  to  change  value  Yp  for  the  value  of  critical  stress  Yc  =  i^Yp  where  cp  <  1  —  coefficient  of  a  longitudinal 
bend  for  centrally  compressed  rod  (projectile)  which  depends  on  its  flexibility  X,  as  well  as  on  properties  of 

Y 

projectile  material.  The  characteristic  type  of  a  dependence  of  normalized  critical  stress  cp  =  —  vs.  flexibility  X 
is  presented  in  Fig. 


Fig.  4  A  characteristic  type  of  a  dependence  of  a  coefficient  of  longitudinal  bend  cp  =  —  vs.  flexibility  X  (the 

Yp 

normalized  diagram  of  a  critical  stress  [30]).  For  values  X>Xe  the  loss  of  stability  occurs  at  elastic 
deformations;  for  values  X<Xe  the  loss  of  stability  occurs  at  elastic-plastic  deformations. 


We  believe,  that  flexibility  of  projectile  can  be  estimated  as  X  =  —  in  the  initial,  most  dangerous  for 

d 

projectile,  moment  of  impact  when  lateral  surface  of  projectile  is  unconstrained  and  the  compressing  stresses  (- 
Go)  appear  in  its  any  cross-section. 

Thus  it  is  easy  to  understand  that  the  loss  of  stability  leads  to  a  decrease  of  value  of  critical  velocity  Usc 
(since  cp  <  1)  (see  (33),  (36)),  especially  for  long  projectiles. 

3.  Numerical  implementation 

We  consider  the  case  of  elastic  deformation  of  the  projectile  and  compare  the  results  with  experimental 
data  where  penetrator  was  only  in  linearly-elastic  state  (it  was  controlled  by  a  post-impact  projectile  state). 

As  an  illustration,  the  penetration  of  projectile  is  calculated  with  a  geometry  given  in  Fig. 5.  Equation 
(1 1)  was  solved  using  Euler’s  method.  For  the  considered  projectile  shapes  and  velocity  range,  the  calculations 
by  formula  (9)  show  that  >  0  on  the  projectile  surface  and  there  is  no  separation  of  material  from  projectile 
(see  Fig.  2).  The  estimate  of  elastic  energy  from  transversal  compression  W\  by  formula  (28)  shows  that  it  can 
be  neglected. 


Fig.  5  Geometry  of  ogival  projectile,  s  =  2>d,  Lh=  1 1.8  mm  for  ogival  head,  s  =  0.5  d,  Li,=  3.555  mm  for 
spherical  head,  Lc  =  59.3  mm,  d  =  7A  \  mm,  L=  Lc  +  Lh. 

In  Fig. 6  the  plots  of  potential  cp  at  the  projectile  tip  changing  with  penetration  depth  are  demonstrated. 
The  dependence  of  derivative  cp^  vs.  P  for  point  J  (Fig.  1)  on  the  projectile  surface  is  shown  in  Fig.7. 


Fig.  6  Dependence  of  potential  at  projectile  tip  on  penetration  depth 
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Fig.  7  Dependence  of  derivative  cpp  vs.  P  for  ogival  (a)  and  spherical  (b)  nose  projectiles 


The  dependencies  of  the  coefficients  A,  B,  C  of  equation  (1 1)  on  penetration  depth  P  are  shown  in  Fig.8. 
Figure  9  shows  the  dependence  of  penetration  velocity  u  vs.  time  t.  Curve  1  here  is  the  solution  of  (1 1)  for  (7o  = 
Us,  curve  2  is  the  solution  of  (11)  with  initial  condition  (26),  curve  3  is  calculated  using  (30').  Figure  10 
demonstrates  the  dependence  of  penetration  velocity  u  vs.  time  t  and  penetration  depth  P  for  ogival  and 
spherical  nose.  The  dependence  of  penetration  depth  P  on  time  t  is  given  in  Fig.ll.  The  typical  fields  of 
velocities  in  target  material  for  spherical  and  ogival  heads  are  illustrated  in  Fig.  12  and  13  for  penetration  depth 
37.0  mm  and  35.9  mm  respectively.  In  the  vicinity  of  projectile  head  these  fields  have  explosive  character.  Rt 
defined  by  (16)  was  used  in  all  calculations  in  Fig.  9 — 13. 
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Fig.  9  Difference  in  dependence  of  penetration  Fig.  10  Dependence  of  penetration 


velocity  u  vs.  time  t  for  ogival  head  at  impact  depth  P  vs.  time  t 

velocity  (7s=  1.0  km/s  for  rigid  (1),  quasi-rigid 

projectile  (2),  and  with  account  for  forced 

elastic  oscillations  (3).  Projectile  —  steel  T- 

200,  target  —  A1  707 1 -T65 1 


t,  |is  P,  cm 


Fig.  11  Dependence  of  penetration  velocity  u  vs.  time  t  and  penetration  depth  P 
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Fig,  12  Velocity  fields  near  the  projectile.  Spherical  head 


Fig.  13  Velocity  fields  near  the  projectile.  Ogival  head. 

Calculations  by  the  model  were  compared  with  available  experimental  data.  Fig.  14  presents  penetration 
depth  dependence  on  impact  velocity  for  projectile  with  spherical  nose.  Target  material  —  aluminium  6061- 
T651,  T=  400  MPa,  E=  68.9  GPa,  v=  0.33,  p  =  2710  IcgW,  projectile  —  steel,  p^=  8000  kg/m^,  Ep  =  210  GPa, 
Vp=  0.29,.  L=  lA.l  mm,  a=  3.55  mm.  The  plots  in  Fig.  15  demonstrate  penetration  depth  dependence  on  impact 
velocity  for  projectiles  with  ogival  nose.  Target  material  in  Fig.  15  (a)  —  aluminium  6061-T651,  Y=  400  MPa, 
E=  68.9  GPa,  v=  0.33,  p  =  2710  kg/m^,  projectile  —  steel  T-200,  Pp=  8020  kg/m^,  Ep  =  210  GPa,  Vp=  0.29,  L= 
82.9  mm,  a=  3.55  mm.  Target  material  in  Fig.  15  (b)  —  aluminium  7075-T651,  Y=  448  MPa,  E=  73.1  GPa,  v= 
0.33,  p  =  2710  kgW,  projectile  —  steel  T-200,  Pp=  8020  kgW,  Ep  =  210  GPa,  Vp=  0.29,  L=  82.91  mm,  a=  3.55 
mm.  Target  resistance  R,  for  calculations  in  Fig.  14,  15  was  evaluated  by  formula  (16).  These  figures  illustrate 
noticeable  influence  of  projectile  forced  elastic  oscillations  on  penetration  depth. 
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Fig.  14  Dependence  of  calculated  penetration  depth  on  impact  velocity  (spherical  nose):  1  -  rigid  projectile,  2  - 

elastic  projectile.  Experimental  points  borrowed  from  [11] 


(a) 


(b) 


Fig.  15  Dependence  of  calculated  penetration  depth  on  impact  velocity  (ogival  nose):  1  -  rigid  projectile,  2 
elastic  projectile.  Experimental  points  were  taken  from:  a)  [7];  b)  [8] 
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Fig.  16  Dependence  of  calculated  penetration  depth  on  impact  velocity  (ogival  nose).  Target  material  — 
limestone,  a  projectile  —  steel  4340.  Experimental  and  initial  data  for  calculation  are  taken  from  [13] 

Fig.  16  shows  comparison  of  model  predictions  with  experimental  data  for  limestone.  Target  resistance 
Rt  was  calculated  by  (19),  p  =  2310  kgW,  projectile  —  steel  4340,  Pp=  7850  kgW,  Ep  =  210  GPa,  Vp=  0.29,. 
L=  71.0  mm,  a=  3.55  mm. 

Results  of  comparison  with  experimental  data  show  reasonably  good  correspondence.  Still  there  is  need 
for  future  investigation  of  oscillation  on  penetration  for  different  pairs  of  projectile-target  both  theoretical  and 
experimental.  The  existence  of  publications  where  good  agreement  of  theoretical  predictions  and  experimental 
data  was  reached  without  involvement  of  oscillation  phenomena  (see  e.g.  Yarin  [16])  may  be  an  indication  that 
a  combination  of  factors  —  from  difference  in  formulations  of  initial  problem  to  improper  account  of  different 
factors  and  their  interaction  —  can  be  the  cause  (and  the  subject  of  further  analysis). 

5.  Conclusion 

New  ballistic  equation  of  normal  friction- free  penetration  in  semi-infinite  targets  of  non-deformable 
axisymmetric  projectiles  of  various  shapes  has  been  proposed  and  analyzed.  It  has  a  simple  and  clear 
mechanical  structure  and  is  built  on  the  base  of  Lagrange-Cauchy  integral  for  non-stationary  irrotational  motion 
of  ideal  fluid.  The  equation  can  be  used  for  analysis  of  penetration  in  elastic-plastic  and  elastic-brittle  materials. 

The  equation  of  projectile  motion  is  constructed  on  the  base  of  definition  of  integral  force  of  resistance 
acting  on  projectile  from  the  part  of  target.  Its  value  is  determined  by  surface  integrals  from  the  functions  of 
velocities  potential  for  target  material  on  the  projectile  surface. 

The  suggested  ballistics  equation  allows  to  make  a  fast  and  rather  reliable  prediction  of  integral 
characteristics  of  normal  penetration  of  projectiles  in  target  (penetration  depth,  drag  force  acting  on  projectile 
and  its  components,  penetration  velocity)  without  additional  fitting  parameters. 

A  number  of  calculations  were  carried  out  for  long  steel  projectiles  of  different  shapes  penetrating 
aluminum  targets  at  velocities  in  the  range  of  600-1600  km/s.  The  results  were  compared  with  available 
experimental  and  computed  data  [7-9,11,13].  A  comparison  of  theoretical  curve  penetration  depth  -  impact 
velocity  with  experimental  data  shows  that  the  experimental  data  are  lower  than  the  theoretical  curve  for  rigid 
projectiles.  The  models  proposed  in  works  [7-9,11,12]  show  a  similar  trend,  and  these  deviations  have  a 
systematic  nature.  In  this  connection  we  took  into  account  forced  elastic  longitudinal  oscillations  of  projectile, 
caused  by  sudden  application  of  braking  force.  The  existence  of  such  oscillations  was  observed  in  recent 
experiments  by  Frew,  Forrestal,  Cargile  [31]  and  also  in  much  earlier  research  by  Masket,  (1949)  [32]  and 
Bluhm  (1956)  mentioned  by  Goldsmith  [28,  chapter  6].  The  energy  of  these  oscillations  caused  by  drag  force 
(work  of  drag  force  spent  on  elastic  oscillations)  is  accumulated  and  has  maximum  in  the  end  of  penetration 
though  deformations  of  penetrator  are  elastic  (small).  It  is  shown  that  these  elastic  oscillations  noticeably  affect 
the  final  depth  and  velocity  of  penetration.  In  our  opinion  longitudinal  oscillations  can  also  reduce  tangential 
stresses  on  the  surface  of  projectile  thus  minimizing  friction  forces.  A  comparison  of  final  penetration  depths 
with  the  experimental  data  is  in  a  good  agreement.  A  critical  impact  velocity  is  defined  at  which  projectile 
undergoes  plastic  deformations  and  loses  its  shape  and  stability  of  motion.  The  obtained  estimate  for  a  critical 
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impact  velocity  has  obvious  mechanical  meaning  and  generalizes  a  well  known  estimate  for  critical  velocity  of 
longitudinal  impact  of  rigid  mass  on  a  rod. 


Appendix 

Example  estimate  of  dynamic  yield  stress 

The  projectile  with  a  spherical  nose  made  from  steel  VAR  4340  has  the  following  parameters 


£'p  =  210GPa,  p^  =  8000  kg/m",  L  =  74.675- 10'^  m,£/  =  7.1  MO”"  m,  J  =  I.S-IO”",  R  =  0.58,  C=  1.00, 
and  the  target  made  of  aluminium  alloy  6061-T65 1  has  the  following  properties 
p  =  2710  kg/m^  Ep  =  68.9  GPa,  v=  0.33,  Y=  400  MPa. 

From  the  experiment  Usc  ~  1 000  ms“’  (see  [9])  and 
2.7M.5  ^ 


,-3 


yp-\^+- 


•74.675 


Tp  =  (1  +  0.0068)  •  Yp  =  1 .0068Tp  GPa, 
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2-400  r 
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689 


=  1455  GPa, 


R;  =1.0-1455  =  1455  GPa,  Cq  = 


210-10"  ,,,,,  _i 

- ^  =  5123  ms  , 

I  8-10^ 


and  from  equation 
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-•5123  +  , 


2-10'’ 


I2.7M0  -0.58 


(l.0068Tp- 1455)  =  1000  ms~\ 


Y  2135 

we  find  Yp  =  2\2>5  GPa.  Quasi-static  yield  stress  is  Yps  =  1481  MPa  (see  [8])  and  ratio  — ^  = - =  1.44,  which 

T„„  1481 

is  a  rather  good  result. 


II  Development  of  new  numerical-analytical  model  of  expansion  of  spherical  cavity  in  brittle 
material  on  the  basis  of  modern  concepts  of  mechanics  of  compressible  porous  and  powder 
materials 


Introduction 

It  is  well  known  that  the  models  of  mechanics  of  expansion  of  spherical  and  cylindrical  cavities  in 
materials  have  essential  significance  for  investigation  of  high-velocity  penetration  in  these  materials.  Both 
models  were  used  for  penetration  problems.  However,  it  is  noted  that  generally  model  of  expansion  of  spherical 
cavity  better  describes  high-velocity  penetration.  Sufficiently  complete  review  and  detailed  analysis  of  existing 
models  of  cavity  expansion  in  different  materials  can  be  found  in  [33].  Notwithstanding  considerable 
achievements  in  construction  of  such  models  it  is  noted  in  [33]  that  sufficiently  simple  and  adequate  analytical 
models  of  cavity  expansion  in  brittle  materials  are  still  in  demand.  Most  widespread  is  the  model  proposed  in 

[23.33] .  In  this  model  for  the  material  in  comminuted  or  pulverized  region  was  used  Mohr-Coulomb  law  (with 
or  without  shear  saturation)  which  defines  the  correspondence  between  compressive  and  shearing  stresses  in  this 
region  using  coefficient  of  internal  friction.  Cavity  expansion  pressure  essentially  depends  upon  this  coefficient 

[23.34] .  This  means  that  material  properties  in  comminuted  region  strongly  influence  the  penetration  resistance 
and  play  important  role  in  its  definition.  To  estimate  this  friction  coefficient  fitting  experiment  is  used  (in 
particular,  normal  impact  of  sphere  on  sample  [34]),  which  essentially  complicates  the  application  of  the  model 
[23,33].  The  new  model  of  spherical  cavity  expansion  in  brittle  material  uses  for  comminuted  region  material 
and  its  boundary  with  cracked  region  rheological  models  of  porous  and  compressible  powder  materials  based  on 
approaches  developed  in  [35-38].  This  allows  to  exclude  friction  coefficient  which  is  determined  in  complex 
experiment.  Instead,  new  parameter,  in  our  opinion,  more  naturally  calculated  parameter  is  introduced:  porosity 
of  fractured  initial  material.  Herein  lies  the  main  difference  of  the  proposed  model  from  existing  ones. 

1.  Main  equations  of  the  model 

Figure  1  presents  geometrical  scheme  of  the  proposed  model  of  expansion  of  spherical  cavity. 
According  to  this  scheme  three  regions  with  different  stress-strain  state  are  formed  in  the  material:  1)  r  >  c  — 
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elastic  region;  2)  b  <  r  <  c  —  dilatation  pore  formation  region;  a  <  r  <  b  —  pulverized  region.  Let's 
consequentially  consider  these  regions. 


Fracture  fronts 


Fig.  1  Scheme  of  the  model  of  cavity  expansion  in  spherical  coordinates  OrcpO: 

1  —  elastic  region  (r  >  c,  volume  deformation  e  =  0);  2  —  dilatation  and  pore  formation  region  (b  <  r  <  c, 
volume  deformation  >  0);  3  —  fractured  material  (a  <  r  <  b,  volume  deformation  e  <  0). 

1.1.  Elastic  region,  r  >  c. 

In  spherical  coordinates  OrcpO  stress-strain  state  in  this  region  is  determined  by  the  following  equations 

[38]: 


MO,.  ^(p  „  ,  . 

— ^  +  2 - =  0 ,  (equilibrity  equation) 

dr  r 

e^-  —  ,  =  60  =  —  (strain-displacement  relations) 

dr  ^  r 

e.  =  ^  [cr,  -  2va^  ]  ,  6^  =  60  =  ^  [( 1  -  v )  -  va,  ]  ,  (Hook's  law) 

with  boundary  conditions: 

m(co)  =  0;  a,p(r  =  c)  =  00(1-  =  c)  =  O/ 


(1) 

(2) 

(3) 

(4) 


where  E,v  —  Y oung's  modulus  and  Poisson  ratio  respectively,  a/ —  material  strength  in  uniaxial  tension. 
Solving  equations  (1) — (4)  we  obtain 


1  +  v  6 

—  ,  a,  =-2a 
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yrj 


-  CT/ 


'  c  ' 


yr  J 


r>c. 


and  for  radial  displacement  u(c)  we  have 


,  X  1  +  v 


(5) 


(6) 


From  (5)  follows  that  in  the  elastic  region  hydrostatic  pressure  p  =  -(a,.  +  a(|,)/3  =  0,  and  volume 
deformation  6  =  6,.  +  26m  =  0. 


1.2.  Dilatation  and  pore  formation  region,  b<r<c. 

Cracked  and  dilated  material  in  this  region  is  formed  at  the  boundary  r  =  c  from  elastic  region  r>  c 
where  tensile  hoop  stresses  Ojp  =  ae  reach  material  strength  at  uniaxial  tension  O/,  i.  e.  boundary  r  =  6  is 
considered  as  a  fracture  front  of  elastic  region  and  front  of  transition  of  material  into  a  state  with  different 
mechanical  properties.  Cracking  and  dilatation  is  accompanied  by  pore  formation.  Thus,  elastic  region  r  >  6  is 
bounded  by  dilatation  and  pore  formation  region  b  <  r  <  c  where  radial  cracks  are  assumed  to  be  distributed  so 
that  hoop  stresses  in  the  whole  region  are  equal  to  zero,  i.  e.  Gjp  =  ag  =  0.  Because  of  this  material  in  the  region  b 
<  r  <  6  is  considered  as  a  continuum  elastically  deformable  in  radial  direction  with  uniaxial  (radial)  stress  state 
with  fracture  occurring  at  the  boundary  r  =  b  where  stress-strain  state  of  material  satisfies  fracture  criterion  and 
material  changes  state  stepwise. 

At  the  boundary  r  =  b  stresses  in  dilatation  pore  formation  region  satisfy  the  following  boundary 
conditions 

Gr(r  =  +  0)  =  a*,  G^{r  =  b  +  Q)  =  G(j{r  =  +  0),  (7) 

where  Gr  —  radial  stresses,  a,,,  ae  —  hoop  stresses,  a  —  stress  defined  below  (see  (21),  (22))  from  fracture 
criterion  at  the  front  r  =  b. 

At  the  boundary  r  =  c  stresses  in  this  region  are  equal 
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aXr  =  c  -  0)  =  -2a/,  a(|,(r  =  c  -  0)  =  a0(r  =  c  -  0).  (8) 

At  crossing  boundary  r  =  c  stresses  ar  are  continuous  and  stresses  acp,  au  are  discontinuous. 

Thus,  stress-strain  state  in  dilatation  and  pore  formation  region  is  defined  by  equations  [9, 1 0] 


^+22^  =  0, 

dr  r 

du 


(9) 

(10) 


with  boundary  conditions  (7)  for  equation  (9)  and  boundary  condition  (6)  for  equation  (10).  Integrating  (9),  (10) 
we  have 


u  = - b 
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'-T 
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(1  +  v) 


(11) 


(12) 


\r  cj 

where  u  =  u(r)  —  radial  displacements  of  material  during  deformation  (u(r)  —  continuous  function  of  argument 

r)- 

From  (8),  (10)  and  (12)  taking  into  account  continuity  of  radial  stresses  at  r  =  c  we  find  position  of 
fracture  front 


c-b, 


(-■) 
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From  (13)  follows  that  dilatation  and  pore  formation  region  exists  for  materials  (c  >  b)  that  satisfy 
condition  2a/ <  (-a  ).  This  condition  is  further  assumed  satisfied. 

From  this  and  from  (12)  for  displacement  u  of  boundary  r  =  b  and  deformations  in  the  considered  region 
at  r  =  +  0  we  have 
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where  e,-  —  radial  deformation;  e^  =  —  hoop  deformations;  Cc  —  volume  deformation,  and  jc  —  intensity  of 

shear  deformations  at  r  =  +  0. 

Because  material  cracking  leads  to  pore  formation  and  volume  deformation  rate  e  is  related  to  porosity 
change  rate  0  by  differential  dependence 
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e  -  - 


or 


de  - 


1-0 

dQ 

1-0 


(15) 


then  material  porosity  0  in  the  dilatation  region  at  r  =  +  0  which  was  formed  in  this  material  after  its  fracturing 

can  be  found  from  the  integral  of  equation  (15) 

1-00 


e  =  ln- 

1-0 

at  e  =  Cc,  00  —  initial  (before  deformation)  porosity  of  material  in  the  region  (b  <r  <  c). 
Thus,  from  ( 1 6)  at  e  =  for  "full"  porosity  0  of  material  at  r  =  +  0  we  have 


(16) 
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0*=l-(l-0o)exp(-e,).  (17) 

From  here  (at  0o  =  0)  for  "induced"  porosity  of  material  0^  (i.  e.  additional  to  0o  porosity  formed  during  material 
cracking  in  the  region  b  <  r  <  c)  we  have 

0c  =  1  -  exp(-ee).  (18) 


1.3.  Pulverized  region  (a  <  r  <  b) 

Region  (a  <  r  <  b)  consists  from  compressed  fractured  material  formed  from  cracked  material  of  region 
(b  <r<  c)  at  the  front  r  =  b  where  material  is  comminuted  (fragmented)  by  shearing  stresses  and  then  densified. 
Because  of  this  we  can  use  for  the  material  of  pulverized  region  rheological  models  of  compactable  porous  and 
powder  materials. 

Further  for  such  materials  we  use  the  following  constitutive  equations  [35-38] 
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where  \|/  =  — 
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,  (p  =  (1  -  0)  —  functions  of  porosity  0;  s  = - index  of  deformation  state  depending  on 
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velocities  of  volume  e  and  shearing  y  deformations;  p  -  — ^  ^ 
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,  Gi,  02,  03  —  main  stresses, 


2  2  2 
(01-02)  +(^2-03)  +(01-03) 


1V2 


;  Fs  —  flow  stress  of  solid  phase  at  uniaxial  compression.  First 


equation  in  (19)  is  a  condition  of  flow  of  porous  body  (in  the  space  of  variables  p  -  x  it  describes  flow  ellipse), 
second  equation  in  (19)  is  an  associated  flow  law. 

At  the  front  r  =  b  (see  Fig.  1)  first  equation  in  (19)  we  consider  as  a  condition  of  fracturing 
(comminution)  by  shearing  stresses  of  a  material  in  the  region  2  with  the  porosity  0  =  0  =1-(1-  0o)exp(-ec), 
see  (17).  Second  —  as  an  associated  law  of  fracturing  (compaction)  of  material  in  the  region  2  with  "induced" 
porosity  0^  =  1  -  exp(-ec),  see  (18)  (i.  e.  porosity  formed  by  radial  cracks  in  region  2  and  additional  to  initial 
porosity  0o).  Flere  7,  is  strength  limit  at  uniaxial  compression  of  porous  free  material.  Thus,  fracturing  and 
densification  of  the  material  in  region  2  at  the  front  r  =  b{i.  e.  formation  of  material  of  region  3)  we  consider  as 
"quasi-plasticity"  of  porous  material. 

Let's  define  the  value  of  deformation  state  index  5  at  the  front  r  =  b.Kt  this  front  it  can  be  written  out  in 
the  form 

'  *  'Ll,  (20) 


y  (/y 
du 


■7 


where  e  =  er  +  2e(p,  y  =  .1—  - - radial  and  tangential  components  of  deformation, 

V  3  '  ^ '  dr  ^  r 

signs  "+",  denote  values  before  front  (at  r  =  b  +  0)  and  behind  the  front  (at  r  =  b  -  0),  respectively.  If  we 
take  into  account  that  <  0,  e<p  >  0,  -e~,  and  =  e  -  2eip,  then  from  (20)  we  obtain 


+  — 
e  -e 
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+  — 
e  -e 


+  + 
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1 

~^(p 

1 

^<p 

+  — 
e  -e 

2  ^  e 


3  e'*'  -e 
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3  e'*'  -  e 
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This  value  of  index  5  =  corresponds  to  compression  in  the  die  mold,  spherical  in  this  case. 

,  0^  +  20  [21  I  f 

If  in  relationships  (19)  we  put  0  =  0,  p- - - — ^ ,  x  =  J—  -  0^| ,  0^  =  0 , 5  =  -J—  (i.  e.  consider 

conditions  (19)  at  the  front  r  =  b)  and  exclude  0<p  from  them,  we  obtain  the  following  equation  for  determination 
of  radial  stress  0.  =  0*  <  0  at  the  front  r  =  b 


0  =-CK, C= 


2(1 -e-) 


3/2 


(21) 


where  porosity  0*  is  calculated  by  formulas  (see  14),(17)). 
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0*  =  l-(l-0o)exp(-ej,  l-(l-v) 

h  1 


(-■) 


It  is  obvious  that  stress  a*  <  0  depends  only  on  material  properties  and  thus  is  its  characteristies. 
Nonlinear  equation  (21)  can  be  solved  by  method  of  consequent  approximations 


-Ys,  0^  =l-(l-0o)exp  l-(l-v) 


K) 


,  ao  =  1(5,  n  =  0,1,2,....  (22) 


It  should  be  noted  that  dimensions  of  regions  2  and  3  (see  (13)  and  Fig.  1)  are  connected  by  relationship 

*  f^<i 

Further  we  assume  that  in  the  region  (a  <  r  <  b)  porosity  0*:  and  volume  deformation  et  are  constant  (do 
not  depend  on  r)  and  also  0^  =  =  0  during  cavity  expansion.  Such  material  state  in  the  region  (a  <  r  <  b) 

takes  place  at  the  fracture  front  r  =  b  (see  Fig.  1).  It  means  that  at  the  level  of  stress  g^=  ct*  material  in  the 
region  (h  <  r  <  c)  at  the  boundary  r  =  b  fractures  by  shearing  stresses,  volume  deformation  and  porosity  of 
fractured  material  in  the  result  of  densification  reach  maximum  values  Ck,  0t «  1 .  After  that  the  material 
becomes  porous  incompressible  medium  which  behaves  in  accordance  with  constitutive  equations  (19)  at 


5  =  (compression  in  spherical  die  mold)  and  0  =  9t.  Because  of  this  from  second  equation  in  (19)  we  have 

[3  0^ 


a,  +  2a„ 


and  intensity  of  shearing  stresses  x  =  ,  —  .  From  this  we  have 

V  3 


G„  1  _  0J-  O'p  ?  O'm  _ 


where  hydrostatic  pressure  p  =  - 


At  the  boundary  r  =  b  takes  place  pressure  and  shearing  stress  discontinuity:  hydrostatic  pressure 

* 

increases  {p{r-b  +  9)-  — ,  p[r-b-9)-a*),  and  intensity  of  shearing  stresses  decreases 


(T(r  =  h  +  0)  =  .^-|a  |,  T(r  =  h-0)  =  y-|G  |0^). 

In  accordance  with  (23)  equilibrium  equation  in  the  region  (a  <r<  b)  takes  form 

— ^  +  2a— =  0,  a-—9i^,aXr  =  b)  =  a, 

dr  r  2 

and  its  solution  is  function 


Thus,  during  cavity  expansion  (i.  e.  with  increasing  radius  a)  region  (a  <  r  <  Z?)  is  joined  at  the  boundary 
r  =  bhy  fractured  and  densified  material  with  porosity  0^. 

Let's  now  estimate  volume  deformation  Ck  and  porosity  0*^  of  the  material  in  the  region  (a  <  r  <  b).  This 
estimate  will  be  done  taking  into  account  elastic  compressibility  of  material  of  solid  phase.  Let's  note  that 
conditions  ( 1 9)  do  not  presuppose  compressibility  of  the  solid  phase,  hence  compressibility  of  the  medium  with 
rheology  (19)  takes  place  due  to  porosity  [35-38]. 

On  the  basis  of  accepted  assumptions,  continuity  equations  in  the  region  (a<r<b)  have  the  form 

Sii  ^ 

— +  2-  =  0  (25) 

or  r 

where  for  convenience  it  is  accepted  t  =  b,  and  from  continuity  of  displacements  at  the  front  r  =  b  for  the 

'¥ 

velocity  li  =  v  of  the  boundary  r  =  b  from  (12)  we  have  expression 
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*  ■/ 1  \  ^ 

V  -u(b)-  —  = - 

^  ^  db 


2a  f 

rJR) 


+■ 


r^^(l  +  v)  -c 


( 2a, 


From  this  we  have  the  following  solution  of  equation  (25) 


u-v 


Krj 


(26) 


(27) 


Let's  define  now  the  relationship  between  average  volume  deformation  ek  in  the  region  (a  <  r  <  b)  and 


ratio  — .  By  definition 
a 


=■ 


8v 


4  3  4  3 

where  5v-—na  —  volume  of  spherical  cavity;  v-—n{^b-u)  —  initial  volume  of  the  region  (a<r<b); 
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u  ,  ^(p  ^ 


1- 


1-v 


2a, 


ih') 


—  elastic  displacement  of  the  boundary  r  =  Zi.  As  a  result  we  have 


relationship 
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(28) 


from  which  it  follows  that  ratio  —  remains  constant  during  cavity  expansion 

b 


1/3 


da 


From  this  and  from  (27)  for  velocity  of  cavity  expansion  =  —  we  obtain 


and 


da 

db 


a 
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(l-4)  i-^k) 
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-  =  const . 


(l-4)  {-^k) 


2/3 


(29) 


at  cavity  expansion  from  zero  radius.  This  relation  and  formula  (28)  define  volume  deformation  Ck  in  the  region 
(a  <r<  b) 


(-;r 


(30) 


Note:  If  we  use  true  deformation,  then  instead  of  formulas  (28)-(30)  we  have,  respectively,  formulas 

1  a  n/3  V* 

- -  =  v  /  , 

b  1-v 


'  a  ' 


1  -  [ajb') 


Porosity  0^  we  find  from  differential  equation 

dQ 


with  initial  condition  e  =  at  0  =  0*.  This  equation  describes  change  of  volume  deformation  e  that  is  caused  by 
change  only  of  porosity  0.  Its  solution  is  function 

0  =  1  -  (1  -  0*)exp(ec  -  e), 

and  for  desired  porosity  we  have  formula 

0  =  1  -  (1  -  0*)exp(ec.  -ek  +  ct^. 


(31) 
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where  —  volume  average  deformation  in  the  region  (a  <r<  b)  caused  by  elastic  deformation  of  solid  phase. 
The  account  of  is  necessary  because  deformation  e*  defined  by  formula  (30)  is  a  sum  of  and  volume 
deformation  =  (ck  -  of  pore  compression. 

The  value  of  can  be  estimated  by  formula 

* 

(32) 

E 

where  — -  —  modulus  of  volume  compression  of  solid  phase  (i.  e.  porousless  material).  Thus, 

3(l-2vJ 

porosity  0*;  is  defined  by  formula 

0  =  1  -  (1  -  0*)exp(ec.  -ek  +  e^),  (33) 


where  e^s  -  — . 


In  the  case  of  small  values  of  0*,  et,  Sh,  formula  (33)  takes  more  simple  form 

^  3(l-v^)  * 

0*  =  0  -ec  +  ek-  ek„  = - - - a  , 


Finally,  from  (24)  for  pressure  Rt  in  the  cavity  we  have 


\aj 


where  a*  is  determined  by  formulas  (21). 


(34) 


Example 

Ceramics  AD995  [1]:  E  =  373.14  GPa,  K  =  231.8  GPa,  a/  =  0.262  GPa,  v  =  0.232,  T,  =  2.26  GPa. 
Equation  (21)  is  solved  by  method  of  simple  iteration  (22): 
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In  the  result  we  have  the  following  figures  for  desired  values: 

a’  =  -10.9  GPa,  =  0.024,  0*  =  0.024,  C  =  4. 15,  v*  =  0,027,  -  =  0.299  ,  «  -0.047,  0^  «  0.018, 

b 


ek  =  -0.029. 


For  the  pressure  in  the  cavity  we  obtain 


7?,  =10.9 


1 


V  0.299  y 


=  11.63  GPa. 


This  value  is  close  to  values  7  —  9  GPa  found  in  literature. 
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III.  Investigation  of  deformation  process  at  high-speed  loadings  of  ceramic  materials  with 
diamond-like  structure  in  adiabatic  approach  at  microlevel 

Introduction 

This  problem  can  be  solved  by  means  of  a  method  of  pseudo-potentials  as  it  was  shown  in  previous  works 
(by  the  example  of  calculation  of  different  physical  values).  In  terms  of  this  method  a  solution  of  the  problem 
can  be  presented  as  follows. 

1.  Calculation  of  pseudo-potential  for  separate  components  of  a  material,  if  a  material  is 
multicomponent. 

2.  Construction  of  a  complete  pseudo-potential  in  view  of  materials  structure. 

3.  Calculation  of  total  energy  of  electron- ion  system  in  terms  of  perturbation  theory  provided  that 
matrix  elements  of  pseudo-potential  of  atoms  are  small  in  comparison  with  Fermi's  energy. 

4.  Stress  appearing  in  a  material  at  any  loading  is  commonly  defined  by  the  first  derivative  of  total 
energy  on  deformation. 

Investigation  of  one  axis  deformation  along  axis  c  in  titanium  by  means  of  interlayer  interaction  showed 
that  densely  packed  basic,  atomic  planes  are  easily  elongated  or  compressed  relatively  each  other.  While 
considering  ideal  structures,  the  basic  planes  remain  parallel  at  one  axis  deformations,  only  distances  among 
those  vary. 

Results  of  investigation  of  process  high-speed  loading  (one  axis)  of  titanium  were  used  as  an  intermediate 
stage  for  transition  to  ceramic  materials  with  diamond-like  structure.  We  described  these  materials  in  hexagonal 
axes  that  enabled  to  allocate  parallel  densely  packed  atomic  planes. 

As  a  result  of  carried  out  numerical  experiment  it  was  obtained  that  densely  packed  atomic  planes  at  high¬ 
speed  one  axis  loadings  are  deformed  in  different  ways  (depending  on  their  distance  from  each  other  and 
relative  positioning  of  atoms  on  these  planes). 

The  developed  technique  is  applied  for  nanoparticles  of  diamond-like  ceramics. 

1.  Calculation  of  mechanical  characteristics  of  diamond-like  materials,  proceeding  from 
energy  of  interaction  of  atomic  planes 

Investigation  of  a  course  of  physical  processes  (deformation,  fracture)  in  materials  between  micro  -  and 
macro  -  scales  remains  an  actual  and  unresolved  problem  of  physical  modeling  in  materials  science  [39]. 

Fracture  process  of  brittle  materials  is  structural  -  multilevel,  and  its  adequate  modeling  is  associated  with 
a  necessity  of  attraction  as  theoretical  methods  and  modern  experimental  and  computer  facilities  [40]. 
Theoretical  methods  which  were  applied  to  a  problem  of  deformation  and  fracture  of  a  wide  class  of  materials, 
originate  from  the  work  of  Barrenblat  [41]  for  quasi-brittle  solid  bodies  (a  model  of  cohesion  zone).  In  terms  of 
continual  theory  the  nonlinear  law  of  boning  of  force  and  displacement  for  atoms  in  top  of  a  crack  is  postulated. 
One  of  the  basic  assumptions  lying  in  a  background  of  all  models  of  cohesion  zone  is  associated  with  a  kind  of 
distribution  of  cohesion  forces. 
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In  terms  of  continual  theory  the  exact  form  of  this  distribution  is  unknown,  but  some  data  can  be  obtained  from 
experimental  works  on  investigation  of  fracture  in  the  field  of  interfaces. 

This  effort  is  to  find  out,  being  based  on  the  method  of  pseudo-potentials,  the  mechanism  of  fracture  at 
microlevel  with  macroscopic  mechanical  characteristics  (diamond-like  materials)  at  one  axis  tension  and 
compression. 

For  a  description  of  fracture  process  of  diamond-like  materials  (diamond  SiC ,  AIN  and  BN )  at  one 
axis  loading  we  use  a  value  of  energy  of  interaction  between  atomic  planes  (perpendicular  to  axis  of  loading).  It 
essentially  differs  from  standard  approach  where  for  calculation  of  mechanical  characteristics  on  microlevel  a 
full  energy  per  atom  in  a  material  is  used  [42-45],  i.e.  on  averaging  on  the  whole  volume,  that  is  not  adequate  to 
a  local  character  of  crack  origin. 

Due  to  anisotropy  of  single  crystal  under  the  influence  of  loadings,  a  material  is  more  easy  deformed  on 
certain  crystallographic  directions  [44-45].  For  investigated  materials  with  cubic  structure  these  directions  are 
body  diagonals  of  elementary  cell.  Atomic  planes,  perpendicular  to  diagonals,  are  densely  packed.  Therefore  it 
is  expedient  to  describe  structure  of  diamond-like  materials  in  hexagonal  axes,  having  chosen  for  an  axis  a 
spatial  diagonal  of  cube  [111].  Thus  the  planes  {111},  being  the  planes  of  octahedron,  would  have  the  indexes 
001  (i.e.  they  are  the  basic  planes). 

Structure  of  diamond  and  diamond-like  materials  would  be  described  in  this  case  by  three-layer 
alternation  of  planes  (001)  like  AA'BB'CC'  where  layers  (planes)  A,  B,  C  consist  of  atoms  of  one  grade 


(  C,  B,  Al,  Si  ),  and  layers  (planes)  A ,  B\  C  -  from  atoms  of  other  grade  (  C,  N)  (see  fig.  1). 


Fig,  1  Description  of  structure  of  diamond-like 
materials  in  hexagonal  axises 

Axis  C  is  represented  by  body  diagonal  of  cube 

[111]  (  c  =  a  Vs  ,  Cl  -  parameter  of  cubic 
lattice).  Atoms  of  one  grade  are  located  on  highly 
densed planes  {111} 
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A  distance  between  two  densely  packed  atomic  layers  A  and  A'  is  d  =  ^  (  C  =  -\/3  ,  - 

parameter  of  cubic  lattice).  Atoms  on  these  layers  are  in  identical  positions.  And  atoms  in  layers  A'  and  B 

(  B'  and  C )  are  displaced  from  each  other  on  Pi  =  —d  +  —b  for  layers  {  A' ,B  ),  or  =  — a +  —b  -  for 

(  B'  ,  O-  Here  a  and  b  -  vectors  of  hexagonal  lattice  on  basic  plane.  A  distance  between  these  layers 
makes  ^2  • 

A  description  of  structure  of  diamond-like  materials  in  hexagonal  axes,  excludes  appearance  of  a 
problem  associated  with  not  identical  deformation  of  materials  on  body  diagonals  [43]. 

Energy  of  interaction  between  atomic  planes  is  calculated  by  means  of  interatomic  potentials  under  a 
scheme  presented  previous  work  [46].  The  aprioristic  pseudo-potential  [47]  is  used  in  calculations.  In  uniform 
model,  from  the  first  principles,  deformations  of  diamond-like  materials  are  presented  at  one  axis  compression 
and  tension  along  direction  [111]. 

The  theory  of  pseudo-potential  can  be  formulated  in  terms  of  pair  potential,  only  in  that  case  when  ions 
are  redistributed  at  constant  average  density  of  a  material. 

By  the  energy  of  interlayer  interaction,  a  theoretical  strength  of  diamond-like  materials  is  calculated  at 
one  axis  tension  and  compression  in  case  of  deformation  [48].  At  small  deformations  when  one  may  neglect  a 
change  of  average  density  of  a  material,  a  calculation  of  theoretical  strength  on  the  basis  of  interatomic  potentials 
is  justified. 

In  the  present  work  we  consider  one  axis  deformation  of  materials  when  a  volume  per  atom  in  a 
material  remains  constant.  Changing  one  of  parameters  of  hexagonal  lattice,  the  second  is  defined  from  a 
condition  of  volume  invariance.  In  this  case  the  employment  of  pair  interatomic  potentials  is  justified  for  any 
possible  deformation. 

A  complete  energy  of  electron  -  ion  system  can  be  written  as 

U=U,+U„  (1) 


where  -  depends  on  volume,  Uy  (energy  of  zone  structure  and  a  part  of  electrostatic  energy)  depends  on 
structure  and  it  can  be  performed  as  a  sum 

=  P) 

energies  of  pair  interatomic  potentials;  i?,  -  interatomic  distance.  Energy  of  zone  structure  defines  indirect 
effect  between  ions 


Ub,  iq{Ri-Rj)  V\q)z{q)e  (q). 


(3) 


Here  E(^)  -  pseudo  potential  of  ions,  ^  -  wave  vector  and  appearing  functions  f(^)and  (  ^  )  are  defined 
in  [9];  i?,  and  Rj  -  radius  vectors  of  ions. 

Energy  of  zone  structure  (3)  together  with  potential  of  direct  interaction  define  a  complete  effective 
interaction  among  ions 
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0(i?)=^  +  -^^fexp(/ (4) 

R  {oTT)  •’ 

where  Z  -  a  number  of  valent  electrons,  0.q  -  atom  volume  in  elementary  cell. 

Let’s  sum  (1.4)  on  two  parallel  highly  densed  layers  distant  from  each  other  on  d  and  shifted  from 
each  other  on  vector  p .  For  the  energy  of  interaction  of  two  layers  (per  square  unit)  we  obtain  [11] 

=  +  — F  exp(z^^J)exp(z^^(.R  +  p))J^.  (5) 

S^(27r)'’l 

Here  S  -  area  per  atom  in  highly  densed  layer;  vector  q  is  divided  on  two  components  -  perpendicular 
highly  densed  plane  and  q^  -  parallel  to  it.  Through  R.  -  we  designate  lattice  vectors  lying  in  a  plane  of  highly 
densed  layer.  To  calculate  integrals  (5)  we  use  equality 

=  (di-gj-  (6) 

^  Si 

Through  -  we  designate  vectors  of  two-dimensional  lattice  inverse  to  considered  plane  of  highly  densed 
lattice  with  constant  a  .  Sum  in  the  right  part  of  equality  (1.6)  equals  zero  if  only  vector  q^  does  not  coincide 
with  one  of  the  vectors  of  reciprocal  lattice. 

If  we  substitute  (6)  in  (5)  we  obtain 

In 

where 

1/  2  7^ 

{gl+qlY^  +_  ,^2 - yT)exp(z^,t/)J^,,  (8) 

L  J  O-YgL+qY 

2n 

and  q^=  —  m;  m=0;  ±1;±2;  ±3 . , 

d 

then 

®(«?,z5)=-|^^^F(J,g^)cos(g^y5).  (9) 

o 

To  calculate  F  {d,g^)  we  have  to  build  the  pseudo-potentials  V  (q). 

Pseudo  potential  of  ions,  accounting  for  materials  structure  in  hexagonal  axises  one  may  represent  as: 

V(q)  =  -Vj+F^  +  Fj  +  F^  +  Fj  +  Fg  ,  (10) 

6 

where  F^,  Fg  -  atomic  pseudo  potentials  of  components  (for  diamond  Fj  =  Fg  ),  5  ^-  ( z  =  1,  5)  radius 
vectors  of  atoms  in  elementary  cell  and  q  -  radius  vector  of  point  of  reciprocal  space.  For  diamond-like 
materials  d  ■  we  obtain  the  following  values: 
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K  -  2^  1-  K  -  2^  1-  7^ 

5  ,=  —c\  5  ^=—a+—h -^—c  ■.  d  ■,=  —a+—b-\ - c; 

4  3  3  3  3  3  12 


(11) 


1 


1 
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d  —a  + —b  + —c  :  d  .=  —a  +  —b^ - c. 

"  3  3  3  '  3  3  12 

After  substitution  of(l.ll)in(1.10),  we  obtain 

r(<?)  =  i 

6 


-_1 

q  c— 

VA+e 

.-,2-  1-  ..1.  ...1.  2-^  ..2. 

— b  -iq—c  -iqi  —  a^ — b)  —iq  —  c 

1+e  ^  ^  e  ^  +e  ^  ^  e  ^ 

_ 

(12) 


and 


V\q)=- 


Va  +VsCOs{q^-) 


x{3  +  2 


c  2  c 

cos  {q^  -  +  A  )  +  cos(- )  +  cos{-q^+q^  {p,  -  p,, ) 


(13) 


Table  1  presents  caleulated  values  of  energy  of  interlayer  interaction  and  lattice  parameters  (in  atomic 
units)  in  equilibrium  state  of  crystals  and  fig.  2  shows  a  dependence  of  stress  vs.  lattice  parameter  cas  a 
material  is  deformed  on  axis  c 

Table  1  Energy  of  interlayer  interaction  O^and  lattice  parameters  ,  Cq  at  equilibrium  state  of  crystal; 
stars  mean  ,  Cq  experimental  values  of  lattice  parameters 


Crystal 

-  Oo- 10  " 

«0 

Co 

* 

* 

Co 

Diamond 

0,25998 

4,7969 

12,007 

4,76 

11,86 

BN 

0,22528 

4,8423 

12,144 

4,89 

12,0 

SiC 

0,23109 

5,735 

14,372 

5,817 

14,23 

AIN 

0,13385 

5,808 

14,549 

5,85 

14,58 

At  known  values  of  energy  of  interaction  of  layers  strain  tensor  is  defined  from  the  ratio  [45] 


cr 


Z 


1  50(<i,y5) 
d  de  ^ 


(14) 


where  -  relative  deformation,  0{d,  p)  -  energy  of  interaction  of  atomic  layers  of  unit  areas  distant  on  d 
and  shifted  relatively  each  other  on  vector  p  . 

Energy  of  interaction  of  atomic  layers  as  A,  A'  -  (  0(^,0))  and  A',  B  -  (  0(^2?P))  differ 

almost  by  an  order.  Theoretical  strength  is  calculated  on  (1.9)  at  energy  0(^2?P)  Uve  times  higher  than 

at  energy  0(^,0).  Thus  in  brittle  materials  with  cubic  structure  and  subjected  to  one  axis  tension  along 

direction  [111],  microcracks  at  first  appear  between  atomic  planes  {111},  taking  place  on  distance  The 

material  most  likely  fails  on  layers  consisting  from  two  atomic  planes  distant  from  each  other  on  ^2  ■ 
Functional  dependence  of  stress  vs.  interplane  distance  for  diamond  is  illustrated  on  fig.  2. 
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Fig.2  Dependence  of  stress  cr  vs.  lattice  parameter  c  for  BN 


A  general  view  of  this  dependence  is  identical  to  the  investigated  materials.  The  interval  BD 
corresponds  to  a  cohesion  zone  (  To  the  left  of  a  point  A  a  function  describes  compression 

process  and  to  the  right  -  tension.  The  point  B  corresponds  to  strength  at  compression,  and  D  at  tension. 

At  very  small  deformations  a  stress-  strain  curve  is  close  to  linear  law,  <j  =  E - .  From  this 

ratio  it  is  possible  to  define  Young’s  modulus  E . 

Values  of  theoretical  strength  <j,  ,  meeting  deformation  at  compression  and  cTj  ,  at 

tension  as  well  as  Young’s  modulus  E  are  presented  in  tab.  2. 

Table  2  Values  of  theoretical  strength  (GPa),  Young’s  modulus  (GPa)  and  strain  rate  at  compression  and 
tension 


Compression 

Tension 

■  Gil 

E, 

max 

^2 

Gil 

^2 

c 

204,4 

0,0975 

1100 

139,07 

0,1456 

1100 

BN 

175,4 

0,0978 

970 

119,24 

0,1458 

970 

SiC 

147,56 

0,0998 

750 

103,12 

0,1459 

750 

AIN 

88,26 

0,098 

450 

59,05 

0,1460 

450 

As  it  can  be  seen  from  the  table,  at  one  axis  tension  or  compression  a  deformation  corresponding  to 
ultimate  strength  of  diamond-like  materials,  possesses  a  certain  value  and  does  not  depend  on  a  material. 
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Surface  energy  of  formed  crack  can  be  calculated  proceeding  from  functional  dependence  of  stress  vs. 
interplane  distance 


r= 


(15) 


c  / 

Here  -  distance  between  basic  planes  without  deformation  (  dQ=  %  ),  and  d^^^  -  corresponds  to 

interplane  distance  when  stress  reaches  its  maximum  value  [12], 

Finally  for  surface  energy  we  obtain  the  following  ratio 

7,78605 -10^  n/m  (16) 

^111 

where  A O  =  O ( d^^^ ,p)-Q>{dQ,p)  -  energy  of  deformation.  The  following  values  of  energy  of 
deformation  are  obtained  at  tension:  A O  =  0,22195  •  10“^  (  0,8081  TO^^)  for  diamond,  0,193595-10^^ 

(  0,698-10”^)  -  for  boron  nitride,  0,189-10^^  (  0,855-10’^  )  -  for  silica  carbide  and  0,113032-10^^ 

(  0,44764- 10~^  )  -  for  aluminium  nitride  (in  brackets  we  put  the  values  of  energy  of  deformation  of 
compression).  =  2^  is  called  as  critical  value  of  rate  of  release  of  deformation  energy. 

With  calculation  thus  there  is  a  transition  from  values  of  microlevel  to  parameters  of  macrolevel. 


Calculation  at  microlevel  enables  to  define  more  precisely  macrocharacteristics  of  materials,  for  example 


Ky^  -  fracture  toughness  coefficient  at  use  of  a  ratio  from  macro-mechanics  [5 1  ] 


Values  Fj^ ,  Ky^  at  tension  and  compression  are  available  in  table  3. 


Table  3  Macroscopic  characteristics  of  diamond 


-like  materials  at  tension  (I)  and  compression  (II) 


Yyc  n/m 

I 

II 

Ky^  MPa  m*'^ 

I  II 

Diamond 

23,6728 

12,839 

5,1067 

3,76075 

3,4-8  [14] 

BN 

20,6484 

11,0914 

4,479 

3,2832 

(  3,5 ±0,5)  [15] 

SiC 

20,1584 

13,5858 

3,89 

3,194 

- 

AIN 

12,056 

7,113 

2,33 

1,79 

- 

Work  [16]  obtain  the  values  Tfj^for  single  crystal  of  cubic  boron  nitride  in  terms  of  (1-4)  MPa*m'^^ 
depending  on  the  amount  of  impurities  content.  The  calculated  values  Fj^  andTfj^,  defined  through 

deformations  and  energy  of  interaction  of  atomic  layers,  are  found  in  terms  of  scattering  of  experimental  values. 
Thus,  the  theoretical  strength  of  diamond-like  materials  in  determined  crystallographic  direction  depends  on  a 
mutual  positioning  of  cooperating  layers.  Densely  packed  atomic  layers  which  are  displaced  relatively  each 
other  (on  vector  Pj  ,  or  p^^ ),  are  more  stable  at  loadings,  than  not  displaced  layers.  The  strength  in  the  first 
case  is  higher  by  an  order  than  in  the  second  one.  Thus,  at  loadings  along  a  direction  [111]  in  diamond-like 
materials  the  cracks  may  appear,  first  of  all,  between  densely  packed  atomic  layers  distant  on  ^ . 

The  calculated  fracture  toughness  coefficient  is  in  terms  of  accuracy  of  experimental  values. 
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Critical  deformation  of  diamond-like  materials  along  direction  [111]  is  a  characteristic  of  structure  and 
does  not  depend  on  a  concrete  material. 

2.  Investigation  of  influence  of  scale  factor  on  theoretical  strength  of  nanoparticle  with 
diamond-like  structure 

A  wide  spectrum  of  practical  applications  of  nanoparticles  as  base  structural  units  of  nanocomposites 
makes  a  theme  of  investigation  of  small  particles  actual.  Since  80*  years  the  interest  to  small  particles  has  brought 
to  a  lot  of  publications  [55-60]. 

A  direct  experimental  investigation  of  physico-mechanical  properties  of  nanoparticles  is  rather  complicated,  as  to 
obtain  authentic  results  is  not  that  simple  methodical  problem  because  of  possible  influence  of  scale  factor  on 
strength  characteristics  of  tested  samples  [56].  Modeling  of  physico  -  chemical  processes  is  complicated  for  these 
systems  by  a  number  of  reasons.  On  the  one  hand,  usual  methods  of  quantum  chemistry  appear  doubtful  in 
application  to  a  description  of  systems  containing  hundreds  of  atoms,  and  on  the  other  hand,  a  macroscopic 
thermodynamics  is  inapplicable  to  small  particles  also  because  of  impossibility  of  division  of  bulk  and  surface 
properties  [57]. 

Among  works  of  the  last  years  a  special  place  is  occupied  by  publications  which  are  devoted  to 
definition  of  surface  energy  depending  on  particles  form  and  size,  with  use  of  elements  of  macro-  and  micro  - 
theories  [57,58]. 

In  work  [58]  to  account  for  surface  energy  of  nanoparticles,  depending  on  their  size  and  forms,  the 
authors  use  a  thermodynamic  model.  Surface  energy  is  expressed  through  a  surface  tension  at  equilibrium 
condition  of  system.  A  consequence  of  given  assumptions  is,  that  all  small  particles  are  compressed  by 
Laplacian  pressure.  Meanwhile,  the  experiment  shows  that  not  always  lattice  parameter  for  crystal  particles 
decreases  with  a  reduction  of  their  sizes.  The  situations  when  lattice  parameter  does  not  vary  even  grows  are 
rather  often  or.  It  means  that  at  reduction  of  linear  sizes  of  a  particle  surface  forces  are  not  always  reduced  to 
Laplacian  pressure.  At  analysis  of  surface  forces  it  is  necessary  to  distinguish  an  equilibrium  case  when  a  body 
has  the  form  corresponding  to  a  minimum  of  system  free  energy  and,  nonequilibrium  when  it  does  not  have  this 
form.  The  body  of  nonequilibrium  form  under  the  action  of  Laplacian  pressure  aspires  to  obtain  the  equilibrium 
form  [57].  In  work  [59]  the  authors  define  surface  energy  of  nanoparticle  through  energy  of  break  of  bond 
between  atoms.  However,  in  this  model  the  calculated  surface  energy  depends  only  weak  on  type  of  structure  of 
nanoparticles. 

In  works  [61-62]  authors  assert  that  a  basic  arbitrariness  in  definition  of  nanoparticle  size  brings  to 
ambiguity  in  definition  of  its  many  mechanical  characteristics. 

A  problem  of  this  investigation  -  is  to  build  (being  based  on  the  method  of  aprioristic  pseudo-potential) 
the  analytical  model  which  adequately  would  describe  mechanical  properties  of  nanoparticles  in  view  of  a  role 
of  advanced  surface  of  particle. 

For  a  description  of  a  process  of  deformation  of  nanoparticles  at  one  axis  loading  in  this  work  we 
suggest  to  use  a  value  of  energy  of  interaction  between  atomic  planes  (perpendicular  to  axis  of  loading). 
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As  nanoparticles  according  to  [63],  possess  properties  of  a  crystal;  thus  let’s  describe  a  structure  of 
nano-diamond  in  hexagonal  axes,  chosen  for  the  axis  a  body  diagonal  of  cube  [111]. 

Energy  of  interaction  between  atomic  planes  is  calculated  by  means  of  interatomic  potentials  which  are 
built  under  above  said  scheme. 

As  models  of  nanoparticle  let’s  use:  a)  nanoplate  with  final  number  of  basic  planes  (  / )  on  spatial  axis 
;  b)  infinite  nanobar  with  a  basis  of  limited  size  on  spatial  axes  X ,  y  ,  and  c)  nanobar  with  limited 
nanosizes  on  three  spatial  axes. 

Earlier  it  was  shown  that  microcracks  first  of  all  appear  between  atomic  planes  {111},  distant  on  ^ ,  in 
massive  brittle  materials  with  cubic  structure  and  subjected  to  one  axis  tension.  The  material  most  likely  fails  on 
layers  consisting  of  two  atomic  planes  distant  from  each  other  on  ■ 

Fig.  3  schematically  shows  a  nanoplate  consisting  of  certain  number  of  atomic  planes  (as  pieces).  Shaded 
strip  correspond  to  strongly  bonded  atomic  planes  distant  on  account  for  those  as  certain  structural 

unit).  Nanoplate  is  represented  as  a  set  of  parallel  structural  units  distant  from  each  other  on  ^ . 

/  =  2 


I 


/  =  4 


I 


/  =  5 


I 


II 


/  =  6 


I 


II 


Fig.  3  Eocation  of  basic  planes  depending  on  their  quantity  / 


At  loading  of  nanoplate,  a  distance  between  neighboring  layers  changes.  In  ideal  case  (infinite  material  on 
three  directions)  the  energy  of  interaction  among  those  is  Og  .  For  materials  of  limited  size  at  calculation  of 

energy  of  interaction  among  layers  it  is  necessary  to  account  for  energy  of  external  surfaces.  In  case  of 
nanoplate  we  have  two  closing  (atomic)  surfaces.  Generally  the  energy  of  interaction  between  neighboring 
structural  units  i.  i  +  \  consists  of  two  parts 

J 


(17) 
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Here  -  a  share  of  energy  of  external  surfaees. 


Assuming,  that  external  surface  of  nanoparticle  possesses  the  energy  which  equals  a  half  of  energy  of 
interlayer  interaction  (with  the  nearest  absent  atomic  plane),  it  is  possible  to  calculate  the  energy  of  interaction 
of  atomic  planes  for  chosen  structural  models. 

If  surface  energy  plays  the  main  role  at  definition  of  mechanical  characteristics  of  materials  thus  there  appears  a 
question  about  a  presence  of  different  values  of  theoretical  strength  of  material  depending  on  distance  up  to  a 
surface.  Closer  to  a  surface  a  material  has  the  higher  strength  than  in  volume.  A  solution  of  this  problem  is 
associated  with  calculation  of  energy  of  interaction  between  structural  units  depending  on  a  distance  of  free 
surface. 

Let's  consider  a  semi-infinite  sample  which  free  surface  coincides  with  highly  densed  plane.  If  for  ideal 
case  (infinite  sample)  the  energy  of  interaction  between  neighboring  structural  units  is  O  ^ ,  because  of  a 
presence  of  energy  of  free  surface  of  a  particle,  its  value  changes.  Let’s  introduce  a  term  "mesolayer",  consisting 
of  two  parallel  structural  units,  a  distance  between  those  makes  ^  4  . 

To  define  the  energy  of  interaction  between  structural  units  at  presence  of  free  surface,  we  suggest  use 
averaging  on  two  neighboring  mesolayers  under  the  following  scheme: 

Let  for  L*  and  2'“^  mesolayer  (fig.  3)  the  energy  of  interaction  in  zero  order  is  Oq  ,  and  a  free  surface  has 
energy  •  Lei’s  define  the  energy  of  interaction  of  L*  meso-layer  as  a  simple  average,  from  ratio 


1  d),  o, 


(18) 


For  2"^^  meso-layer  the  energy  of  interaction  we  define  falling  from  refined  value  of  energy  of  L'  meso-layer 


and  for  j  layer 


(19) 


O  ,  = 


v+l 


(20) 


Thus  a  distribution  of  surface  energy  on  meso-layers  is  obtained  depending  on  distance  of  external  surface. 
If  the  object  has  limited  size  on  spatial  axis  z  then  it  is  necessary  to  account  for  the  influence  of  energy  of  2"'* 

O, 

surfaces.  Let  the  second  surface  has  the  energy  .  The  method  of  calculation  of  distribution  of  energy  of 
interaction  in  view  of  external  surface  is  the  same,  as  above  described.  As  a  result  we  obtain  for  i  layer 


o  +0),  — 4 

!■  0  1  2  ry  j+2 


7-2 


(21) 


where  j  -  a  number  of  layers  in  nanoparticle.  A  number  of  layers  j  and  a  number  of  atomic  layers  /  are 
connected  with  each  other 


36 

1  =  2  j  ,  then  Oj  =  Oj, 

1  =  2  j  + 1 ,  respectively,  ®2“®o- 

The  last  item  in  (2.5)  provides  the  law  of  preservation  of  full  energy  of  system  at  limited  number  of  atomic 
planes  /  .  The  ratio  (2.5)  can  be  used  for  calculation  of  energy  of  interlayer  interaction  both  for  bulk  samples 
and  for  nanoplates  as  well. 

In  case  of  infinite  nanobar  (the  basic  plane  of  nanoparticle  has  sizes  n  a  and  m  b  (fig.  4),  and  the 
amount  of  those  is  infinite,  n,  m  -  integers),  we  account  for  energy  of  lateral  surfaces. 


Fig.4.  Basic  plane  of  diamond-like  materials  in  hexagonal  axises  a,  b ,  axis  c  is  perpendicular  to  the 

plane 


According  to  problem  condition,  we  calculate  mechanical  characteristics  at  one  axis  loadings,  using  the 
energy  of  interaction  between  structural  units.  The  last  are  perpendicular  to  the  axis  of  loading.  Therefore,  in  lieu 
of  lateral  surfaces  we  consider  lateral  atomic  lines  on  planes.  The  energy  of  interaction  between  lateral  atoms, 
located  on  neighboring  planes,  represents  a  required  share  of  surface  energy.  Thus,  the  energy  of  interaction 
between  atomic  layers  (2.1)  can  be  presented  as 

0  =  (iV,Oo+iV2®.)/(A^i+A^2),  (22) 


where  -  a  number  of  atoms  on  planes  located  inside  a  nanobar,  N2  -  extreme  atoms  of  these  planes,  O  ^  - 
energy  of  interaction  between  extreme  atoms  in  neighboring  layers.  In  this  case 

=  (n  - 2) (m  - 2)  and  N2=2{m  +  n-2)  (23) 

Accounting  for  a  number  of  bonds  of  extreme  and  angle  atoms  of  nano-diamond  on  basic  planes  for  energy  of 
meso-layer  interaction  we  obtain 

0=0o  l  +  A2  +  f)-^  ,  (24) 

3  m  n  9m  n 


when  a  basic  plane  is  a  parallelogram  with  sides  MU  and  n  b  ,  and 
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^  ^  1  1/1  ^  5 

O  — Oq  Ih — ( — I — )3 - 

3  m  n  3m  n 


for  a  rectangle  with  sides  fnci  3L-  and  n  b  and  with  equal  basic  square. 

2 

If  m,  n ^  00,  then  O  ^  Oo . 

In  case  of  limited  bar  or  nanobar  for  i  layer  we  obtain 


^  1  1 .  1  lx  8 

Here  0^=0^  ^  +  +  - 

3  m  n  9mn 


^  ^  1  1/1  lx  5 

O  — Og  1h — ( — I — )h - 

3  m  n  3m  n 


(^=0,  1,2) 


At  known  values  of  energy  of  meso-layer  interaction,  a  strain  tensor  is  defined  from  the  ratio 

..4^,  (29, 

a  oe  ^ 

where  e  ^  -  relative  strain  rate,  O  (  J  )  -  energy  of  interaction  of  atomic  layers  in  elementary  cell  distant  on  d  . 
Dependence  of  theoretical  strength  on  distance  up  to  external  surface  of  nano-diamond  (  i  -a  number  pf  meso- 
layer)  from  size  of  a  basic  plane  (  m,  n),  from  nanobar  height,  or  nanoplate  thickness  (  /  ,  at  fixed  number  of 


meso-layer)  is  illustrated  in  fig.  5-7. 
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Fig.5  Dependence  of  theoretical  strength  of  diamond  nanoparticle  vs.  distance  from  external  surface  (vs.  a 

number  of  mesolayer  i ,  /  =  1 20 ;  m  =  n  =  25 ) 
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Fig.  6  Dependence  of  theoretical  strength  of  diamond  nanoparticle  vs.  size  of  basic  area  (vs.  m,  n),  when 

/  =  120,  a  i  =  3 


a  I  - Qj 


40 


Fig.  7  Dependence  of  theoretical  strength  of  nano-diamond  vs.  length  (or  /  -  a  number  of  atomic  layers) 
nanobar,  when  I  -  m  =  n  =  250 ,  and  II  -  m  =  n  =  25  An  both  cases  i  =  3 


In  a  case  of  nanobar,  consisting  from  120  highly  dense  packed  atomic  planes,  on  distance  from  a  surface 
*  fz  50  t 

=  a  v3  —  (  a  -parameter  of  cubic  lattice),  value  of  theoretical  strength  at  compression  and  tension  is 

6 


minimum  (see  fig.  3).  At  even  values  /  a  theoretical  strength  varies  depending  on  distance  up  to  external 
surface  as  symmetric  function  that  can  not  be  said  about  nanoparticle,  containing  odd  number  of  densely  packed 
atomic  layers  (see  fig.  1).  At  odd  number  of  atomic  layers  external  surfaces  of  nanoparticles  possess  different 
energy  that  results  in  asymmetrical  distribution  of  energy  on  nanobar  length.  At  increase  of  a  number  of  atomic 
planes  a  common  picture  of  distribution  of  energy  does  not  vary.  In  case  /  =  1200,  on  distance  from  the  first 


surface 


h  =a*  Vs  —  and  from  the  second  surface  h-,  =a*  Vs — 


hi ,  a  value  of  theoretical  strength  does  not 


vary  and  coincides  with  a  value  of  theoretical  strength  for  massive  samples  (tab.  4). 

Table  4  presents  values  of  theoretical  strength  at  compression  cTj  and  tension  cTj  for  diamond-like 
nanoparticles  as  nanobar,  depending  on  its  size  and  form  of  basis  (rectangular  or  parallelogram). 


Table  4  Value  of  theoretical  strength  (GPa)  for  diamond-like  nanoparticles  depending  on  their  size  and  form 
of  its  basis  (basic  plane).  Here  /  -  a  number  of  atomic  layers,  m,n  -  integers  responsible  for  size  of 
nanobar  basis  and  i  -  a  number  of  meso-layer  describing  a  distance  to  external  surfaces 


Nanobar  with 

rectangle  basis 

Nanobar  or  nanoplate  with  parallelogram  basis 

1  =  100;  m  =  n  =  25 

i  =  2 

/  =  100;»j  =  n;  1  -  120;  m  -  n  -  2500 

i  =  2  30<f<570 

0-2 

0-2 

0-2 

c 

236,709 

161,048 

235,768 

160,4 

204,4 

139,07 

BN 

203,14 

138,087 

202,3 

137,5 

175,4 

119,24 

SiC 

170,88 

119,42 

170,203 

118,94 

147,56 

103,12 

AIN 

102,208 

63,385 

101,802 

68,113 

88,26 

59,05 

Let's  note,  that  at  values  m  =  n  =  2500 ,  theoretical  strength  at  compression  and  tension  for  nanoplate  with  the 
basis  of  rectangular  or  parallelogram,  do  not  differ  from  each  other.  Investigation  (computing  experiment) 
showed  that  at  values  m  =  n  >  500  a  value  of  theoretical  strength  (at  tension  and  compression)  does  not  depend 
on  the  form  of  basic  area.  Comparative  calculations  by  an  example  of  SiC  nanobar  with  the  basis  of  rectangular 
and  parallelogram  is  presented  in  tab.  5. 


Table  5  Dependence  of  theoretical  strength  of  SiC  nanobar  vs.  the  form  of  basic  area  (I-  parallelogram,  II 
rectangle) 


Form 

m  =  n 

m  =  n 

m  =  n 

m  =  n 

m  =  n 

m  =  n 

m  =  n 

=  25 

=  50 

=  100 

=  150 

=  200 

=  300 

=  500 

I 

0-1 

180,023 

179,845 

178,709 

178,323 

178,129 

177,934 

Ml, 111 

I 

0-2 

127,207 

125,68 

124,89 

124,62 

124,486 

124,349 

124,24 
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II 

0-1 

182,749 

180,027 

178,75 

178,34 

178,14 

177,939 

177,777 

II 

^2 

127,715 

125,80 

124,92 

124,63 

124,49 

124,353 

124,24 

3.  CONCLUSION 


-  Strength  of  nanoparticles  at  one-axis  deformations  strongly  depends  on  particle  size  along  the  axis  of 
loading  and  distances  from  external  surface  of  particle; 

-  Theoretical  strength  of  nanobar  (nanoplate),  at  even  value  of  atomic  planes,  along  the  axis  of 
deformation  on  identical  distances  from  external  surfaces,  has  identical  values.  It  is  associated  with  identical 
energy  of  free  surfaces  of  nanoparticle,  unlike  nanoparticles  with  odd  number  of  planes.  For  nanoplate  with 


size 


h  >aS—, 

3 


a  value  of  theoretical  strength  on  distance 


—  from  external  surface  always  aspires  to  a 


value  for  massive  samples. 

-  Theoretical  strength  does  not  depend  on  the  form  of  nanobar  basis  of  at  values  of  basis  area 

5  >44200  a 


IV.  Ballistic  limit  velocity.  Theory  and  Experimental  Approximation  on  the  Base  of  Indentation 
Technique  (A4  New  Empirical  Law  for  Critical  Velocity  of  Penetration) 


TABLE  OF  SYMBOLS 

Pp  -  density  of  KEP  material; 

-  density  of  target  material; 

Lp  -projeetile  length; 

u  -  veloeity  of  movement  of  a  eontaet  surfaee; 

R  -  statie  penetration  resistanee  of  target  material  (defined  by  analytieal  solution  of  a  problem  of 

expansion  of  spherieal  eavity); 

t  -  time,  s; 

Ug  -  impaet  veloeity; 

-  ballistie  limit  veloeity; 

hf  -  target  thiekness; 

Y  -  yield  stress  of  target  material; 

HV  -  Viekers  hardness  of  target  material; 

HM  -  Meyer  hardness  of  materials; 

djj  -  plastieity  eharaeteristie  of  target; 

dp  -projeetile  diameter; 

nip  -projeetile  mass; 

bottom  indexes  t  and  p  respeetively  are  referred  to  target  and  projeetile. 
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INTRODUCTION 

In  a  problem  of  interaction  of  kinetic  energy  projectile  (KEP)  with  a  target  a  calculation  of  ballistic  limit 
velocity  represents  a  significant  scientific  and  theoretical  interest.  A  process  of  penetration  of  indenter  into 

target  has  much  in  common  with  a  process  of  indentation  of  a  material  with  rigid  indenter.  In  this  work,  being 
based  on  the  representations  about  indentation  process  elaborated  by  authors  the  problem  of  calculation  of 
ballistic  limit  velocity  is  submitted  on  the  basis  of  experimental  results  obtained  at  indentation  (microhardness 
HV  and  plasticity  characteristic  ).  The  theory  of  calculation  V^j.  for  a  case  of  not  deformable  projectile  is 
offered  as  well.  For  the  common  case  (deformable  and  not  deformable  projectile)  empirical  equations  which 
enable  to  calculate  with  satisfactory  accuracy  are  suggested 
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1.  Dynamic  tests  of  materials 

In  this  work  we  investigated  a  wide  class  of  target  materials:  aluminum  alloys  (standard  and  high- 
strength,  developed  in  IPMS  NAS  of  Ukraine),  titanium  alloy,  steels  (soft  and  strengthened)  and  ceramics.  A 
choice  of  objects  of  investigations  is  stipulated  by  a  necessity  to  obtain  the  materials  with  a  wide  range  of 
physico-mechanical  properties.  It,  in  its  turn,  was  necessary  to  obtain  rather  representative  experimental  data  file 
on  the  basis  of  which  it  would  be  possible  with  a  rather  good  reliability  to  define  regularities  for  ballistic  limit 
velocity  of  penetration.  Thickness  of  targets,  geometrical  sizes  and  a  material  of  projectile  were  varied  as  well. 

A  chemical  composition  of  the  investigated  alloys  is  presented  in  table  1,  and  physicomechanical 
properties  -  in  table  2. 

Samples  for  tests  represented  the  plates  of  100x100  mm  of  various  thicknesses,  which  were  fixed  in 
holder’s  slot.  The  holder  was  steel  thermo-treated  plate  (HRC  40)  with  thickness  of  40  mm.  The  plate  has 
square  reach-through  hole  of  50x50  mm  for  free  flight  of  kinetic  projectile  in  case  of  penetration  of  tested 
sample.  The  holder,  in  its  turn,  was  rigidly  installed  in  the  frame  of  machine  (Fig.l).  With  the  help  of  this  holder 
the  sample  can  be  located  in  required  position  i.e.  to  provide  a  hit  of  projectile  in  certain  place  of  the  sample 
(with  an  accuracy  ±  5mm)  and  under  a  certain  angle  to  sample  plane  (in  this  work  this  angle  was  90°). 


Fig  1  Equipment  for  dynamic  tests 


A  distance  from  a  channel  cut  of  missile  machine  to  a  sample  plane  made  10  m.  Projectile  velocity  was 
fixed  by  means  of  serially  produced  device  "Polet"  with  an  accuracy  of  0.1  m/s.  A  principle  of  velocity 
registration  -  fixing  of  time  of  projectile  flight  between  two  photo  sensors  established  on  a  distance  of  0.5  m. 

With  the  purpose  to  increase  a  reliability  of  test  results  a  thick  polyethylene  film  was  mounted  behind  a 
tested  plate  on  distance  of  30  mm.  The  plane  of  film  tension  was  in  parallel  to  a  plane  of  plate.  The  fact  of  plate 
penetration  was  considered  as  defined  at  presence  of  reach-through  hole  in  the  film. 
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2.  Theoretical  estimation  of  ballistic  limit  velocity 


Calculation  of  ballistic  limit  velocity  of  a  material  for  a  case  of  not  deformable  projectile  (when  its 

hardness  considerably  exceeds  the  hardness  of  target)  is  suggested.  A  theoretical  substantiation  is  given  below. 

The  equation  of  projectile  movement  in  infinite  target  is  considered  and  it  looks  as: 

du  1  2  (  \ 

PpLp  —  =  --PtU  -R,  u[Q)=u,.  (1) 

dt  2 

This  equation,  in  essence,  is  a  special  case  of  the  equation  of  V.P.  Alekseevskii  [1]  if  in  the  last  to  assume  that 
projectile  is  not  deformable. 

Equation  (1)  is  an  equation  with  divided  variables  and  its  solution  can  be  presented  as 


Pt 


1 


-t  = 
ILpPp  a 


Me  M 

arctg  —  -  arctg  — 

V  a  a  J 


2 

a  =  — 
Pt 


or 


Us  -  a -tg 


u{t)=  a 


a  -  ugtg 


f  \ 

apd 

X^^pPp  j 
f  \  ’ 

apd 

\^RpPp  j 


a  = 


2R 

Pt 


For  penetration  depth  P  from  equation  (1)  we  obtain 


P  = 


t  u 

j* u{T)dT  “  “ 


'^PpRp  du  _  PpRp 

2  2  2  2 
U  +  M  Pf  U  +  M 


(2) 


(3) 


If  from  (2),  (3)  to  exclude  M  ,  then  we  obtain  a  dependence  of  penetration  depth  vs.  time  t 


p(t)= 


p  pRp 

Pt 


+ul\a-tgyf 
i^^a-tgyf+{us-tgy)^ 


7  = 


aPd 

'^RpPp 


Ballistic  limit  velocity  is  approximately  determined  from  (3)  as  impact  velocity  m^  (=  ) ,  at  which 

P  =  hf  and  m  =  0  (i.e.  projectile  stops),  where  -  target  thickness.  Therefore  from  (3)  we  obtain 

PpRr 


ht  = 


d^ln 

Pt 


''  ^ 
2 


V 


a 


(4) 


from  this  equation  for  velocity  F^^  we  have 


=  a  exp 


(  u  ^ 

pA 

-1  = 

2R 

exp 

[  pA  ^ 

-1 

ypp^p  J 

Pt 

'\Pp^p ) 

.(5) 
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Introducing  a  designation  a  =  — ^  let’s  re-write  formula  (5)  as: 

Pp^p 


= 

^  cr 


2R 

Pt 


(5’) 


For  ideal  elastic-plastic  material  for  R  one  may  use  formula  [22] 

r  =  y[ 


2  ,  IE 

-I-  In 


(6) 


,3  ZY ) 

At  the  same  time  R  may  be  estimated  as  the  function  of  HV .  In  work  [3]  it  was  shown,  that  at 
estimation  of  target  resistance  against  penetration  of  a  kinetic  energy  projectile  it  is  necessary  to  account  for 
plasticity  characteristic  5]-[  determined  (as  hardness  HV),  at  indentation  of  a  material  by  Vickers’s  indenter 


[4],  A  product  of  resistance  of  target  material  R  on  plasticity  characteristic  5fj  (or  product  HM  ■  )  enables 

to  allocate  a  part  of  resistance  associated  with  plastic  deformation  of  a  material  at  penetration  of  kinetic  energy 
projectile. 

In  a  view  of  the  above  said,  3  variants  of  estimation  R  are  carried  out  with  the  aim  of  achievement  of 
best  agreement  of  theory  and  experiments  (table  2). 
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Table  1  Chemical  composition  of  target  material 


Material 

groups 

Material 

markiug 

Chemical  composition,  mass.% 

HV, 

GPa 

s„ 

Al-alloy 

6112 

Al-lMg-0.3Cu-0.6Si 

0.84 

0.92 

AMr-5 

A1  -5Mg 

0.80 

0.93 

8112 

Al-8Zn-2.lMg-l.9Cu-0.15Zr 

1.80 

0.84 

8312 

Al-8Zn-2.1Mg-  1.9Cu-0.15Zr-0.12  Sc 

1.9 

0.82 

Steel 

Steel  124 

Fe  -  0.38C  -  0.25Si  -  0.035Cr  -  1.6Mn  -  0.013  Ni  -  0.41Mo 

3.7 

0.85 

Steel  600 

Fe  -  0.68C  -  0.29Si  -  0.61Cr  -  0.77Mn  -  1.87  Ni  -  0.32Mo 

6.72 

0.71 

Steel  40 

Fe  -  0.38C  -  0.36Si  -  0.08Cr  -  0.9Mn  -  0.12  Ni  -  0.12Mo  - 
0.07CU 

2.02 

0.922 

High  hardness 
steel 

Fe  -  0.25C  -  0.36Si  -  0.27Cr  -  0.9Mn  -  1  .ONi  -  0.3Mo  - 
0.1 9Cu 

3.45 

0.866 

Ti-alloy 

Ti-alloy 

Ti  -  5.1V  -  0.3Fe  -  0.4Si  -  3.0A1 

3.22 

0.854 

Ceramics 

SiC 

Self-bounded  silica  carbide 

22 

0.35 

Metallic  plates 
(dams), 
used 
for 

experimeuts 
with  ceramics 

HI 

Steel  45  raw  (TP  free) 

1.81 

0.93 

n2 

Steel  45  (hardening  -I-  tempering) 

3.22 

0.88 

n3 

AMG5 

0.80 

0.93 

n4 

High  hardness  steel 

Fe  -  0.25C  -  0.36Si  -  0.27Cr  -  0.9Mn  -  1  .ONi  -  0.3Mo  - 
0.1 9Cu 

3.45 

0.866 
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Table  2  Results  of  dynamic  tests  of  materials  for  the  case  of  «not  deformable  projectile» 


Material 

class 

Material 

Target 
thickness 
h, ,  mm 

dp, 

mm 

d-‘p  5 

Mm 

Pt^ 

g/sm^ 

HV, 

GPa 

^Ht 

HVp 

GPa 

^Hp 

<^0.2  ’ 
GPa 

A(3) 

GPa 

Y  actual 
^  cr 

Calculation  by  formula  (5) 

ycalcul 

m/s 

Vcfcul^2) 

m/s 

ycalcul 

m/s 

Al-alloys 

6112 

8 

7.62 

26 

2.7 

0.84 

0.92 

1.6 

0.93 

0.26 

1.406 

251 

326 

340 

AMG5 

5 

9 

12 

2.7 

0.80 

0.93 

1.6 

0.93 

0.25 

1.372 

387 

290 

379 

393 

8112 

8 

7.62 

20.5* 

2.8 

1.8 

0.84 

1.9 

0.93 

0.56 

2.503 

438 

418 

522 

569 

8312 

8 

7.62 

20.5* 

2.8 

1.9 

0.82 

1.9 

0.93 

0.59 

2.561 

456 

425 

524 

579 

AMG5 

5 

7.62 

16 

2.7 

0.80 

0.93 

1.78 

0.93 

0.25 

1.372 

290 

249 

325 

337 

Steels 

Steel  40 

40.5 

12.7 

51.6* 

7.8 

2.02 

0.922 

6.55 

0.746 

0.4 

2.388 

845 

784 

854 

889 

High 

hardness  steel 

16.5 

\1.1 

51.6* 

7.8 

3.45 

0.866 

6.55 

0.746 

0.95 

4.660 

685 

558 

671 

721 

Steel  45 

3 

7.62 

12.5* 

7.8 

1.81 

0.93 

1.6 

0.93 

0.5 

2.409 

470 

355 

409 

424 

Ti-alloy 

Ti-alloy 

17 

12.7 

51.6* 

4.6 

3.22 

0.854 

6.55 

0.746 

1.0 

2.849 

490 

526 

515 

557 

*-  for  given  materials  in  calculations  we  used  core  length  but  not  projectile  as  jacket  fails  at  penetration 


R{\)  =  HM-dH 

'2 


r{i)  =  y\ 

R{2>)=Y 


+  In 


+  In 


3YJ 

2E^ 

vsrj 


■Sh  {y  -  (Jq  2  ) 


a 


PT^t 
P  pL  p 
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The  best  coincidence  of  theory  and  experiment  gives  the  variant  (2)  in  which  the  parameter  is  accounted 
for.  In  this  case  experimental  and  theoretical  results  (obtained  without  fitting  parameters)  almost  coincide. 

As  it’s  already  been  said  above,  the  suggested  theoretical  estimation  is  applicable  only  for  a  case  when 
projectile  hardness  considerably  exceeds  hardness  of  target.  It  is  well  seen  at  the  analysis  of  results  for  alloys 
8112  and  8312  where  projectile  and  target  hardness  are  close.  Divergences  of  theory  and  experiment  here  are 
more  essential. 

3.  Empirical  estimation  of  ballistic  limit  velocity 

For  more  common  case  of  an  estimation  of  ballistic  limit  velocity  (including  a  deformable  projectile) 
semi-empirical  formulas  are  suggested. 

On  the  basis  of  available  in  literature  and  obtained  by  authors  results  the  following  empirical  equations 
for  parameter  are  suggested: 


C=C,A- 

(7) 

P 

vl=C,p,'^(HVS„), 

(8) 

P 

(9) 

P  (HVSh)p 

Cl,  C2  5  ^*3  are  the  empirical  constants  determined  by  available  experimental  results  (in  our  case  -  according 

to  tests  of  10  different  materials)  by  determination  of  average  arithmetic  value.  P  is  so-called  static  pressure  of 
projectile  on  target. 

m„ 

P  = - 

/  4 

where  and  dp  -  projectile’s  mass  and  diameter  respectively. 

Experimental  results  and  calculations  under  formulas  (7),  (8)  and  (9)  are  presented  in  table  3. 

Formula  (7)  is  most  simple  and  demands  a  minimum  of  initial  data  for  calculation  which  are  easy  to 

obtain  by  direct  measurement  and  weighing.  However,  a  coincidence  of  calculated  and  experimental  data  is  not 
enough.  The  error  reaches  40  %  and  higher  (for  steel  600  in  particular).  Though  for  ceramics  and  high-strength 
aluminum  alloys  (8112  and  8312)  a  quite  good  conformity  is  observed. 

Formula  (8)  provides  as  a  whole  the  higher  accuracy  (especially  for  aluminum  alloys)  though  the  values 
of  micromechanical  properties  of  target  (hardness  and  plasticity  characteristic)  are  necessary  for  calculation.  For 
ceramics  this  formula  does  not  work  (an  error  up  to  80%). 

And  at  last  there  comes  formula  (9).  For  calculation,  besides  micromechanical  properties  of  target, 
knowledge  of  projectile  properties  is  required.  Calculation  under  this  formula  gives  the  best  agreement  of 
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calculated  and  actual  values  .  The  error  of  calculation  does  not  exceed  5-6%  for  metals  and  8-10  %  for 
ceramics  and  only  in  some  cases  reaches  15-20%  for  ceramics  SiC  (the  explanation  of  this  is  given  below). 

The  analysis  of  the  above-stated  empirical  dependences  enables  to  draw  the  following  conclusion.  For 
the  most  correct  calculation  it  is  necessary  to  account  for  mechanical  properties  of  both  target  and  projectile 

(besides  their  geometrical  sizes  as  well  as  their  density  and  mass).  Therefore  further  it  is  expedient  to  use 
formula  (9)  for  calculations,  as  it  provides  the  greatest  accuracy. 

It  is  should  be  noted  that  at  calculation  for  composite  targets  "ceramics  -  dam"  a  ratio  was  used: 


=V^  +  Vu 

^  comp  ^  cer  '  ^  bai 


where  and  ballistic  limit  velocities  of  penetration  of  ceramic  target  and  dam 


respectively. 
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Table  3  Dynamic  test  results 


Materials 

Class 

Alloy 

A 

g/cm^ 

ht 

mm 

nip 

g 

dp 

mm 

mm 

HV, 

GPa 

HVp 

GPa 

^Hp 

y  actual 
^  cr 

m/s 

ycalcul 

empirical 

m/s 

yCalcul  ^2^ 

empirical 

m/s 

ycalcul^^ 

empirical 

m/s 

6112 

2.7 

8 

7.9 

7.62 

26 

0.84 

0.92 

1.6 

0.93 

323 

435 

276 

330 

AMTS 

2.7 

5 

5.9 

9 

12 

0.8 

0.93 

1.6 

0.93 

387 

474 

295 

353 

8112 

2.8 

8 

3.6* 

7.62 

20.5* 

1.8 

0.84 

1.9 

0.927 

438 

457 

405 

447 

8312 

2.8 

8 

3.6* 

7.62 

20.5* 

1.9 

0.82 

1.9 

0.927 

456 

457 

411 

453 

< 

AMTS 

2.7 

5 

5.5 

7.62 

16 

0.8 

0.93 

1.80 

0.93 

290 

412 

256 

289 

Steel  124 

7.8 

4.5 

4.8* 

7.62 

26* 

3.7 

0.85 

1.78 

0.931 

689 

508 

649 

737 

Steel600 

7.8 

4.5 

4.8* 

7.62 

26* 

6.72 

0.71 

1.78 

0.931 

916 

508 

799 

906 

SteeMO 

7.8 

40.5 

30.1* 

12.7 

51.6* 

2.02 

0.922 

6.55 

0.746 

845 

1219 

1200 

794 

High  hardness  steel 

7.8 

16.5 

30.1* 

12.7 

51.6* 

3.45 

0.866 

6.55 

0.746 

685 

779 

970 

642 

n  1 

7.8 

3 

4.8* 

7.62 

26* 

1.81 

0.93 

6.55 

0.746 

<280 

412 

385 

255 

n  1 

7.8 

3 

2.3* 

7.62 

12.5* 

1.81 

0.93 

1.80 

0.93 

470 

542 

507 

574 

n2 

7.8 

3 

2.3* 

7.62 

12.5* 

3.22 

0.88 

6.55 

0.746 

328 

412 

500 

331 

n4 

7.8 

4.5 

3.6* 

7.62 

20.5* 

3.45 

0.866 

6.55 

0.746 

- 

582 

726 

480 

n4 

7.8 

4.5 

30.1* 

12.7 

51.6* 

3.45 

0.866 

6.55 

0.746 

- 

406 

507 

335 

Ti-alloy 

Ti-alloy 

4.6 

17 

30.1* 

12.7 

51.6* 

3.22 

0.854 

6.55 

0.746 

490 

606 

725 

480 

SiC 

3.09 

8 

4.8* 

7.62 

26* 

22 

0.35 

6.55 

0.746 

- 

423 

846 

560 

SiC 

3.09 

6 

3.6* 

7.62 

20.5* 

22 

0.35 

6.55 

0.746 

- 

423 

846 

560 

SiC 

3.09 

10 

30.1* 

12.7 

51.6* 

22 

0.35 

6.55 

0.746 

- 

381 

763 

505 

y 

SiC 

3.09 

12 

30.1* 

12.7 

51.6* 

22 

0.35 

6.55 

0.746 

- 

418 

836 

553 

s 

SiC 

3.09 

15 

30.1* 

12.7 

51.6* 

22 

0.35 

6.55 

0.746 

- 

467 

934 

618 

r_) 

SiC  +  n  1 

8+3 

4.8* 

7.62 

26* 

6.55 

0.746 

660 

590 

930 

615 

SiC  +  n  2 

8+3 

4.8* 

7.62 

26* 

6.55 

0.746 

710 

590 

983 

650 

SiC  +  n  3 

8+5 

4.8* 

7.62 

26* 

6.55 

0.746 

660 

526 

869 

606 

SiC+n4 

6+4.5 

3.6* 

7.62 

20.5* 

6.55 

0.746 

730 

719 

1115 

738 

*  -  for  given  materials  in  calculations  we  used  length  and  mass  of  core  but  not  projectile,  as  jacket  fails  at  penetration; 
Projectile’s  density  Pp  in  all  cases  was  accepted  as  7.8  g/cm^. 

Values  of  empirical  constants  Cj ,  C2  ^nd  C3  were  1.52x10^  m^/s^,  0.79x10^  m/kg  and  1.69x10^  m^/s^  respectively. 
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This  equation  we  obtained  from  the  law  of  conservation  of  energy.  It  was  supposed,  that  kinetic  energy 
of  projectile  was  consistently  spent  on  penetration  through  ceramic  target  and  dam.  Projectile’s  mass  in  a 
process  of  penetration  was  considered  as  constant.  It,  probably,  not  absolutely  meets  a  reality.  Projectile  while 
going  through  ceramics  loses  a  part  of  mass  and  only  the  rest  part  of  projectile  influences  on  a  dam.  Therefore  a 
real  value  )  must  be  higher  than  calculated  one).  It  takes  place. 

Let’s  get  back  to  formulas  (7)  -  (9).  Approximately  while  considering  a  projectile  as  a  cylinder,  let’s 
mip 

find  its  volume  F  = - L  ,  where  L  is  the  projectile’s  length.  Density  of  projectile’s  materials  is: 


m. 


m. 


Pp  = 


V  Td," 


Since 


m. 


nd". 


=  P,  then  Pp=P  Lp  or  P  =  PpLp.  By  introduction  of  designation  a ,  used  in 


theoretical  formula  (5’): 

pA 


a  = 


Pv^ 


p  p 


formulas  (7)  -  (9)  one  may  write  as: 

vl=qa 

rl  =  C2a(Hrs„), 

(HVSh)^ 

The  empirical  dependences  written  down  in  this  form  obtain  some  physical  sense.  Let’s  decipher 
formula  (9  ’),  as  the  most  correct. 

Let’s  note,  that  parameters  a  in  formulas  (7’)  -  (9’)  and  (5)  are  completely  similar.  And  accounting  for 
that  at  a  «  1  (that  takes  place  in  our  case)  exp  a  —  1  =  a  a  certain  similarity  between  the  above  said  formulas 
is  observed.  A  difference  of  formula  (9’)  is  that  it  accounts  for  projectile  properties  that  are  required  for  our  case 
when  a  projectile  can  be  deformed. 

Also  it  should  be  noted  that  the  authors  of  this  project  had  introduced  a  concept  of  material  resistance  to 
penetration,  determined  as  HV  ■  Sfj  [64].  Thus  (HV  ■  5^^  is  the  resistance  against  penetration  of  target 


(7’) 

(8’) 

(9’) 


material  (which  is  used  here  in  lieu  of  that  one  calculated  by  formula  (6)).  By  analogy  let’s  introduce  a  concept 
(HV^Sh)  as  the  penetrating  ability  of  projectile  material. 

In  view  of  it,  it  is  possible  to  say  that  a  square  of  ballistic  limit  velocity  is  directly  proportional  to 
specific  surface  mass  of  target  and  resistance  to  penetration  of  target  material  and  inversely  to  specific  surface 
mass  and  penetrating  ability  of  projectile  material  (formula  (9’)). 


52 


The  constant  in  equation  (9’)  is  numerically  equal  to  a  square  of  ballistic  limit  velocity  of  penetration  of 
a  target  by  projectile  made  from  the  same  material  and  with  a  length  equal  to  target  thickness. 

4.  Comparative  characteristic  of  dynamic  resistance  of  different  materials 

Being  based  on  the  obtained  experimental  results  and  formula  (9)  the  estimation  of  dynamic  resistance 
of  investigated  materials  is  carried  out.  As  a  criterion  a  critical  value  ( )cr  of  target  material  (representing  its 
specific  surface  mass)  at  which  the  material  withstands  an  influence  of  projectile  of  certain  type  moving  with  a 
speed  Vp  is  accepted.  Along  with  it  a  critical  thickness  of  target  material  is  defined.  Results  of  calculations 

are  presented  in  table  4. 

The  obtained  data  witness  that  the  least  dynamic  resistance  is  attributed  to  soft  aluminum  alloys.  Further 
come  high-strength  aluminum  alloys  and  soft  steel  40.  Then  (in  terms  of  increase  of  dynamic  resistance)  we 
have  titanium  alloy  and  high-strength  steel.  Then  -  super  high-strength  steel  600.  And  at  last  ceramics  SiC 
which  possesses  the  best  dynamic  resistance  from  the  investigated  materials.  It  should  be  noted,  however,  that  in 
practice  the  advantage  of  ceramics  is  not  that  impressive.  The  reason  is  that  ceramics  requires  so-called  "dam" 
(usually  -  a  steel  plate),  providing  for  a  required  rigidity  for  tested  sample  and  interfering  its  fracture 
immediately  from  a  dynamic  impact.  A  presence  of  dam  brings  a  parameter  to  increase  for  a  composite 

target  "ceramic  +  dam"  (see  table  4). 
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Table  4  Dynamic  resistance  of  materials  against  penetration  of  different  projectiles 


Target 

material 

Projectile 

material 

HV, 

GPa 

^Ht 

Pp 

g/sm^ 

Lp 

mm 

Pt 

g/sm^ 

Vp  =700  m/s 

Vp  =1000  m/s 

Vp  =1500  m/s 

{Ph)tcr 

g/sm^ 

cr 

mm 

{ph)tcr 

g/sm^ 

cr 

mm 

{Ph)tcr 

g/sm^ 

6112 

O) 

o 

II 

a 

S 

c^ 

*  ^ 
o 

II 

a 

> 

a 

0.84 

0.92 

7.8 

26 

1.67 

2.7 

12.7 

47.2 

26.0 

96.3 

58.5 

216.6 

AMG5 

0.8 

0.93 

7.8 

26 

1.67 

2.7 

13.2 

49.0 

27.0 

100.0 

60.8 

225.0 

8112 

1.8 

0.84 

7.8 

26 

1.67 

2.8 

6.5 

23.3 

13.3 

47.4 

29.9 

106.8 

8312 

1.9 

0.82 

7.8 

26 

1.67 

2.8 

6.3 

22.6 

12.9 

46.0 

29.0 

103.6 

Steel  124 

3.7 

0.85 

7.8 

26 

1.67 

7.8 

3.1 

4.0 

6.4 

8.2 

14.4 

18.4 

Steel  600 

6.72 

0.71 

7.8 

26 

1.67 

7.8 

2.1 

2.6 

4.2 

5.4 

9.5 

12.1 

Steel  40 

2.02 

0.92 

7.8 

26 

1.67 

7.8 

5.3 

7.6 

10.8 

13.8 

24.3 

31.1 

High  hardness 
steel 

3.45 

0.87 

7.8 

26 

1.67 

7.8 

3.3 

4.2 

6.7 

8.6 

15.1 

19.4 

Ti-alloy 

3.23 

0.85 

7.8 

26 

1.67 

4.6 

3.6 

7.7 

13 

15.8 

16.4 

35.6 

SiC 

22 

0.35 

7.8 

26 

1.67 

3.1 

1.3 

4.2 

2.6 

8.5 

5.9 

19.1 

6112 

Hardened  steel, 

(HVp  =  6,55GPa,  5hp=0.75 

0.84 

0.92 

7.8 

26 

4.88 

2.7 

37.1 

137.5 

75.8 

280.6 

170.5 

631.2 

AMG5 

0.8 

0.93 

7.8 

26 

4.88 

2.7 

38.6 

142.8 

78.7 

291.5 

177.1 

655.9 

8112 

1.8 

0.84 

7.8 

26 

4.88 

2.8 

19.0 

67.8 

38.7 

138.3 

87.1 

311.2 

8312 

1.9 

0.82 

7.8 

26 

4.88 

2.8 

18.4 

65.8 

37.6 

134.2 

84.6 

302.0 

Steel  124 

3.7 

0.85 

7.8 

26 

4.88 

7.8 

9.1 

11.7 

18.6 

23.9 

41.9 

53.7 

Steel  600 

6.72 

0.71 

7.8 

26 

4.88 

7.8 

6.0 

7.7 

12.3 

15.7 

27.6 

35.4 

Steel  40 

2.02 

0.92 

7.8 

26 

4.88 

7.8 

15.4 

19.8 

31.4 

40.3 

70.7 

90.7 

High  hardness 
steel 

3.45 

0.87 

7.8 

26 

4.88 

7.8 

9.6 

12.3 

19.6 

25.1 

44.1 

56.5 

Ti-alloy 

3.23 

0.85 

7.8 

26 

4.88 

4.6 

10.4 

22.6 

21.3 

46.2 

47.8 

103.8 

SiC 

22 

0.35 

7.8 

26 

4.88 

3.1 

3.7  (4.8*) 

12.1 

7.6  (9.8*) 

24.7 

17.1 

55.6 

6112 

o 

II 

& 

■■  £ 

^  to 

G  O 

G  VO 

S  II 

a 

> 

a 

0.84 

0.92 

7.8 

51.6 

4.88 

2.7 

73.7 

272.9 

150.4 

557.0 

338.4 

1253.2 

AMG5 

0.8 

0.93 

7.8 

51.6 

4.88 

2.7 

76.5 

283.5 

156.2 

578.5 

351.5 

1301.7 

8112 

1.8 

0.84 

7.8 

51.6 

4.88 

2.8 

37.7 

134.5 

76.9 

274.5 

172.9 

617.7 

8312 

1.9 

0.82 

7.8 

51.6 

4.88 

2.8 

36.6 

130.5 

74.6 

266.4 

167.8 

599.4 

Steel  124 

3.7 

0.85 

7.8 

51.6 

4.88 

7.8 

18.1 

23.2 

36.9 

47.4 

83.1 

106.6 

Steel  600 

6.72 

0.71 

7.8 

51.6 

4.88 

7.8 

11.9 

15.3 

24.3 

31.2 

54.8 

70.3 

Steel  40 

2.02 

0.92 

7.8 

51.6 

4.88 

7.8 

30.6 

39.2 

62.4 

80.0 

140.4 

180.0 

High  hardness 
steel 

3.45 

0.87 

7.8 

51.6 

4.88 

7.8 

19.1 

24.4 

38.9 

49.9 

87.5 

112.2 

Ti-alloy 

3.23 

0.85 

7.8 

51.6 

4.88 

4.6 

20.6 

44.9 

42.1 

91.6 

94.8 

206.1 

SiC 

22 

0.35 

7.8 

51.6 

4.88 

3.1 

7.4  (9.6*) 

24.0 

15.1  (19.6*) 

49.0 

34.0 

110.3 

*  The  values  for  the  eomposite  target  with  the  ceramies/back  thickness  ratio  3:1  are  shown  in  parenthesis 
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5.  Maps  of  distribution  of  microhardness,  plastic  deformation  and  plasticity 
characteristic  in  targets  after  the  impact  loading 

For  investigation  of  deformation  mechanism  and  strain  hardening  during  impact  loading  and 
penetration,  targets  from  alloys  with  different  types  of  strengthening  were  used:  steel  40,  high  hardness 
steel  and  Ti-alloy.  Chemical  compositions  of  alloys  are  given  in  the  table  1. 

Conditions  of  impact  loading  of  targets  from  alloys  are  given  in  the  part  1 . 

For  investigation  of  the  plastic  deformation  around  the  penetration  channel  the  incompletely 
penetrated  targets  were  selected.  The  targets  were  cut  along  the  axis  of  moving  of  the  kinetic  energy 
projectile  (KEP),  and  the  Vickers  microhardness  (at  a  load  of  2  N)  was  measured  in  the  cut  section. 
Typical  maps  of  microhardness  distribution  are  given  in  Fig.2,  4  and  5. 

A  modeling  of  the  deformation  process  was  carried  out  by  compression  of  specimens,  which  were 
prepared  from  the  same  steels  and  Ti-alloy.  The  specimen  size  was  5x5x5  mm^.  The  deformation  curves 
in  the  coordinates  compression  stress  (Jcomp  ~  pl^stic  deformation  Sp  are  given  in  Fig.6.  The 

microhardness  was  measured  in  the  specimens,  which  were  deformed  to  different  deformation  degrees. 
The  dependence  of  microhardness  on  the  deformation  degree  is  given  in  Fig.  7  as  well. 

The  local  plasticity  characteristic  of  the  material  5h  [65-67]  was  calculated  for  every  point,  in 
which  microhardness  was  determined. 

This  plasticity  characteristic  has  been  introduced  as  a  part  of  the  plastic  deformation  in  the  total 
elastic-plastic  deformation  under  the  indenter 

=  (10) 

Flere  Sp ,  and  are  the  mean  values  of  plastic,  elastic  and  total  deformation  of  the  specimen  on  the 
contact  area  of  the  indenter  with  specimen  in  the  loading  direction;  =  Sp  +  s^  and  Sp  ~  const  when 
a  pyramidal  indenter  is  used  for  hardness  measurement. 

As  it  can  be  seen  from  Eq.(lO),  5fj  is  a  dimensionless  parameter  ranging  from  0  (for  “pure” 
elastic  deformation)  to  1  when  the  elastic  deformation  is  insignificant.  It  was  shown  in  [68],  that  elastic 
deformation  is  given  by 

(11) 

E, 

where  HM  is  the  Meier  hardness,  Ei  is  the  Young  modulus  and  Kj  is  the  Poisson’s  ratio  of  the 
material. 

Plastic  deformation  can  be  calculated  from  the  ratio  [66-68]: 


55 


Sp  =-  ln(sin  f  2  A  (12) 

where  /2  the  angle  between  the  face  of  the  hardness  indentation  pyramid  and  the  loading  direction.  In 
this  case  ^2  >  7i  >  where  is  an  appropriate  angle  of  diamond  indenter.  The  angle  T2  1°^  pyramidal 


indenters  can  be  defined  from  the  ratio  [66,  69] 


ctgy2  =ctgyi  -1.77- 


Here  is  the  effective  Young’s  modulus  of  the  indenter  -  specimen  contact  couple,  and 

1  \-v}  1-v^ 


where  E2  and  V2  are  Young’s  modulus  and  Poisson  ratio  of  the  indenter. 

The  values  of  plasticity  characteristic  for  investigated  materials  are  given  in  the  Figs. 2,  4  and 
5.  Using  the  dependence  E[V(spj  (Fig.7)  as  a  calibration  curve,  the  value  of  the  conditional  deformation 
Sp  was  calculated  for  every  point  of  target  as  well.  Consequently,  the  maps  of  microhardness  distribution 
in  Fig.2,  4  and  5  are  simultaneously  the  maps  of  plastic  deformation  distribution  and  distribution  of  the 
plasticity  characteristic  S fj .  The  distribution  of  microhardness  along  the  lines  Oa,  Ob  and  Oc  of  target 


from  steel  40  is  given  in  Fig.3,  the  lines  {Oa,  Ob  and  Oc)  are  plotted  in  Fig.2. 
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Area 

HV ,  GPa 

8  p ,  % 

5h 

2.8  ^3.4 

20  +  >40 

0.89  +  0.86 

2.5  ^2.8 

8  +  20 

0.90  +  0.89 

1 - 1 

1.2 -r  15 

1.5  +  8 

0.91+0.90 

1 - 1 

2.0 +  2.2 

0.5+  1.5 

0.92  +  0.91 

1 _ 1 

1.9 +  2.0 

0  +  0.5 

-0.92 

Fig  2  The  target  from  steel  40,  velocity  of  KEP  V  =  420  m/s. 

The  map  of  distribution  of  microhardness  HV  ,  the  value  of  plastic  deformation  £p 

and  plasticity  characteristic  after  impact  loading. 
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Area 

HV ,  GPa 

Sp ,  % 

Sh 

4.5  ^4.7 

>60 

0.83+0.82 

1  1 

3.8  +4.5 

15  +  60 

0.86  +  0.83 

1 - 1 

3.6 +  3.8 

1.5  +  15 

0.86  +  0.86 

3.2 +  3.6 

0+1.5 

0.88  +  0.86 

Fig  4  The  target  from  high  hardness  steel,  veloeity  of  KEP  V  =  468  m/s. 

The  map  of  distribution  of  mierohardness  HV  ,  the  value  of  plastie  deformation  £p 


and  plasticity  characteristic  after  impact  loading. 
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Area 

HV ,  GPa 

Sp ,  % 

Sh 

3.8  +4.2 

7  +  22 

0.81+0.79 

I  I 

3.3  +3.8 

1-rl 

0.84  +  0.81 

- 1 _ I - 

3.2 +  3.3 

0  +  2 

-0.84 

Fig  5  The  target  from  Ti-alloy,  velocity  of  KEP  V  =  274  m/s. 

The  map  of  distribution  of  microhardness  HV  ,  the  value  of  plastic  deformation  £p 


and  plasticity  characteristic  after  impact  loading. 


0  10  20  30  40  50  60 

8  ,% 

P 

Fig  6  The  influence  of  deformation  degree  Sp  on  the  compression  stress  (Jcompp 
for  static  compression  loading  for  different  materials  of  targets. 
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0  10  20  30  40  50  60 

^  ,  % 

P 

Fig  7  The  influence  of  deformation  degree  Sp  on  the  microhardness  HV 
for  static  compression  loading  for  different  materials  of  targets. 

The  value  of  plastie  deformation  during  impaet  loading  of  a  target  by  the  KEP  ean  be  estimated 
from  the  Fig.2,  4  and  5.  It  is  seen  that  plastie  deformation  is  essential,  and  may  be  equivalent  to  40-50  % 
of  statie  plastie  deformation  for  steels. 

In  the  ease  steel  40  the  deformed  zone  has  more  dimension  in  eomparable  with  zones  for  high 
hardness  and  Ti-alloy.  The  hardness  of  the  deformed  zone  inereases  and  aehieves  the  value  HV  =3.4 
GPa  for  steel  40  (Fig.2). 

As  it  is  seen  from  Fig.2,  4  the  inereasing  of  hardness  around  penetration  ehannel  for  investigated 
steels  aehieves  50  %,  that  eorresponds  to  plastie  deformation  £'^=40-60  %.  For  target  from  Ti-alloy  the 

plastie  deformation  around  penetration  ehannel  is  lower  (£^»20  %)  and  at  this  deformation  fraeture  is 
observed  in  dynamie  and  in  statie  loading  as  well. 

6.  The  energy  expended  for  the  plastic  deformation 

The  method  for  estimation  of  the  energy  Wp  expended  for  the  plastie  deformation  of  target 

during  impaet  loading  or  penetration  of  KEP  was  proposed  by  authors  in  [64]. 

The  energy  expended  for  the  plastie  deformation  of  the  target  during  penetration  of  the  indenter 
may  be  estimated  for  every  i-region  in  the  map  of  Fig.2  as 

Wi  =  iTiRiSiUi . 


(15) 
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Here  is  the  surface  of  i-region  in  Fig.2  with  approximately  constant  hardness  is  the 

distance  from  the  axis  of  the  penetration  channel  to  the  centre  of  gravity  of  the  i-region;  is  the  energy 
expended  for  plastic  deformation  of  a  unity  of  volume  in  the  i-region. 

It  was  assumed  that  Uj  may  be  evaluated  by  equation 

^  p  i  -  - 

Ui  =  J  CTidSp  =  ,  (16) 

0 


where  and  Spi  are  the  average  values  of  hardness  and  plastic  deformation  for  the  i-region. 

The  total  energy  expended  for  plastic  deformation  of  target  will  be 

.  (17) 


w  = 
P 


Calculation  of  w„  were  fulfilled  for  steel  40.  The  kinetic  energy  of  KEP  was  calculated  as  well 


mv 

- )• 

2 


Table  5  The  energy  expended  for  the  plastic  deformation  of  targets 
during  penetration  of  the  indenter  and  kinetic  energy  of  the  KEP. 


Material 

V  ,  m/s 

WpJ 

Wp  /w^,% 

Steel  40 

420 

2655 

1338 

~50 

It  is  seen  that  the  value  of  the  energy  expended  for  the  plastic  deformation  of  target  during 
penetration  of  the  KEP  Wp  is  equal  approximately  50  %  of  the  kinetic  energy  of  the  KEP  w*.  Other  50  % 
of  the  kinetic  energy  is  expended  for  the  heating,  elastic  vibration  etc. 

7.  CONCLUSION 

1.  The  theory  enabling  to  calculate  ballistic  limit  velocity  of  penetration  of  target  by  not 
deformable  projectile  is  developed.  Calculated  and  experimental  data  are  in  a  good  agreement. 

2.  The  semi-empirical  formula  allowing  to  calculate  ballistic  limit  velocity  both  for  not 
deformable  and  deformable  projectile,  at  different  ratio  of  mechanical  properties  of  projectile  and 
target  is  suggested.  Values  of  hardness  and  plasticity  characteristic  8h  defined  at  indentation  and 
introduced  by  authors  are  employed  in  calculations  for  characteristic  of  mechanical  properties.  As 
targets  we  investigated  different  metals  (steels,  aluminum  and  titanium  alloys)  as  well  as  ceramics. 
As  projectile  materials  we  applied  soft  and  hardened  steels. 

The  obtained  experimental  results  are  in  a  good  agreement  with  calculations  made  under  the 
suggested  formula. 

The  estimation  of  critical  thickness  of  target  for  different  types  of  projectile  and  various  velocities 
of  their  movement  is  carried  out. 
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3.  By  critical  specific  surface  weight  of  targets  the  investigated  materials  can  be  arranged  in  the 

following  order  in  line  with  growth  of  density  of  target  per  surface  unit  (when 

=  const ): 

•  soft  aluminum  alloys, 

•  high-strength  aluminum  alloys 

•  soft  steel, 

•  titanium  alloy, 

•  high-strength  steels, 

•  ceramics 

4.  The  maps  of  distribution  of  microhardness,  plastic  deformation  and  plasticity  characteristic  in 
targets  from  aluminum  alloys  and  steels  after  impact  loading  were  obtained. 

The  energy  expended  for  the  plastic  deformation  around  the  KEP  in  the  case  of  its  stopping  in  the 
steel  targets  (F  <  )  is  approximately  50  %  from  the  KEP  kinetic  energy. 
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